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Abstract 

The Chiral Quark-Soliton Model of nucleons is based on two ideas: 1) the major role 
of spontaneous chiral symmetry breaking in hadron physics and 2) the relevance of the 
large N c (= number of colours) limit for the real world. In these lectures I review the 
theoretical foundations of the model, the physics involved, and some of applications. 
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1 How do we know chiral symmetry is spontaneously 
broken? 



The QCD lagrangian with Nf massless flavours is known to posses a large global symmetry, 
namely a symmetry under U(Nf) x U(Nf) independent rotations of left- and right-handed 
quark fields. This symmetry is called chiralf\. Instead of rotating separately the 2-component 
Weyl spinors corresponding to left- and right-handed components of quark fields, one can 
make independent vector and axial U(Nf) rotations of the full 4-component Dirac spinors - 
the QCD lagrangian is invariant under these transformations too. 

Meanwhile, axial transformations mix states with different P-parities. Therefore, were 
that symmetry exact, one would observe parity degeneracy of all states with otherwise the 
same quantum numbers. In reality the splittings between states with the same quantum 
numbers but opposite parities are huge. For example, the splitting between the vector p and 
the axial a\ meson is (1260 — 770) ~ 500 MeV; the splitting between the nucleon and its 
parity partner is even larger: (1535 — 940) ~ 600 MeV. 

The splittings are too large to be explained by the small bare or current quark masses 
which break the chiral symmetry from the beginning. Indeed, the current masses of light 
quarks are: m u ~ 4 MeV, rrid ~ 7 MeV, m s ~ 150 MeV. The only conclusion one can 
draw from these numbers is that the chiral symmetry of the QCD lagrangian is broken down 
spontaneously, and very strongly. Consequently, one should have light (pseudo) Goldstone 
pseudoscalar hadrons - their role is played by pions which indeed are by far the lightest 
hadrons. 

The order parameter associated with chiral symmetry breaking is the so-called chiral or 
quark condensate: 

(-0-0) ~ —(250 MeV) 3 . (1.1) 

It should be noted that this quantity is well defined only for massless quarks, otherwise it is 
somewhat ambigious. By definition, this is the quark Green function taken at one point; in 
momentum space it is a closed quark loop. If the quark propagator has only the 'slash' term, 
the trace over the spinor indices implied in this loop would give an identical zero. Therefore, 
chiral symmetry breaking implies that a massless (or nearly massless) quark develops a non- 
zero dynamical mass (i.e. a 'non-slash' term in the propagator). There are no reasons for 
this quantity to be a constant independent of the momentum; moreover, we understand that 
it should anyhow vanish at large momentum. The value of the dynamical mass at small 
virtuality can be estimated as one half of the p meson mass or one third of the nucleon mass, 
that is about 

M(0) ~ 350 - 400 MeV; (1.2) 

this scale is also related to chiral symmetry breaking and should emerge together with the 
condensate (|1.1|). 

One could imagine a world without confinement but with chiral symmetry breaking: it 
would not be drastically different from what we meet in reality. There would be a tightly 

2 The word was coined by Lord Kelvin in 1894 to describe moleculas not superimposable on its mirror 
image. 
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bound light Goldstone pion, and relatively loosely bound p meson and nucleon with approx- 
imately correct masses, which, however, would be possible to 'ionize' from time to time. 
Probably the spectrum of the highly excited hadrons would be wrong, though even that is 
not so clear |]J. We see, thus, that the spontaneous chiral symmetry breaking is the main 
dynamical happening in QCD, which determines the face of the strong interactions world. 
If one understands the microscopic mechanism of spontaneous chiral symmetry breaking 



and knows how to get the quantities (|1.1|J1.2| ) from the only dimensional parameter there is 



in massless QCD, namely Aqcd> one gets to the heart of hadron physics. My sense is that it 
is achieved by ways of the QCD instanton vacuum, see || |3|] for recent reviews. Therefore, 
I start by showing in section 2 how instantons lead to a low-energy theory exhibiting chiral 
symmetry breaking and the appearance of a momentum-dependent constituent quark mass 
M(k). In section 3 the resulting effective chiral lagrangian is theoretically studied and in 
section 4 it is applied to build the Chiral Quark-Soliton Model. 



2 Low-energy limit of QCD from instantons 

The idea that the QCD partition function is dominated by instanton fluctuations of the gluon 
field, with quantum oscillations about them, has successfully confronted the majority of facts 
we know about the hadronic world (for a review see ref. 0). Instantons have been reliably 
identified in lattice simulations (for a review see ref.[|]]), and their relevance to hadronic 
observables clearly demonstrated || (]. Thus, the effective low-energy theory coming from 
instantons seems to be well-motivated. The aim of this section is to derive this effective 
theory to which QCD is reduced at low momenta. 



2.1 Some of the results 

There is a well-known general statement that to get chiral symmetry breaking one needs a 
finite spectral density i>(\) of the quark Dirac operator at zero eigenvalues, since the chiral 
condensate is proportional to exactly this quantity: (ipip} = — nv(0)/V^ @. A natural way 
to get z/(0) ^ is to have a finite density of instantons and antiinstantons (J's and J's for 
short) in the 4-dimensional space-time, N/V {4) . Indeed in the presence of the topologically 
non-trivial gluon fluctuations fermions necessarily have an exact zero mode P, as it follows 
from the Atiah-Singer index theorem. In the ensemble of J's and J's the would-be zero 
modes in the background field of individual J's and J's are smeared into a band with a finite 
spectral density at zero eigenvalues flOfl , leading to z/(0) ^ 0. The instanton vacuum provides 



thus a beautiful mechanism of chiral symmetry breaking JTT| . 

There are two mathematically equivalent ways to treat quarks in the instanton vacuum. 
One is to calculate an observable in a given instanton backgound and then to average over 
the collective coordinates of J's and J's and sum over their total numbers, N+ and AL. This 



approach has been developed in refs. [11], [T|| . The quark propagator in the instanton vacuum 

3 It can be added that in the solvable N = 2 supersymmetric version of QCD it is instantons — and 
nothing besides them — that seem to be sufficient to reproduce the expansion of the exact Seiberg-Witten 
prepotential 0. 
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takes the form of a massive propagator with a dynamically generated momentum-dependent 
mass (in Euclidean space, hence the factor i in the 'non-slash' term): 



W = £^S- ^Hconst.^J^-), (2.1) 

Here N/V is the instanton density at equilibrium and p is the average instanton size, F(pp) 
is a combination of modified Bessel functions and is related to the Fourier transform of the 
would-be zero fermion mode of individual instantons, 



F(pp) = 2z 



/„(z)tfi(z) - h(z)K (z) - -h{z)K x {z) 



- VP 12 P~*°° (pp) 3 



° F(0) = 1. (2.2) 



The numerical constant in eq. ( J2.1| ) is of the order of unity and is determined by the self- 
consistency or gap equation: 

C J (27T) 4 M2(p)+p2 V 

The chiral condensate (ipip) is the quark propagator taken at one point; in momentum 
space it is given by a quark loop: 



- WW = «*V>« « 4iV c / ^JjS^p = —t' • (2.4) 

To get the numerical estimates of the condensate and of the constituent quark mass one 
may rely on the variational calculation of the instanton vacuum characteristics ||, |13[ , which 

(3) 

relates them to the only dimensional parameter in QCD, Kqcd- Taking = 280 MeV, 
one finds from refs. |9|, [13j the basic characteristics of the instanton vacuum, namely the 
average distance between neighbouring instantons, R = (N/V)" 1 ^ and their average mean 
square radius, p, to be 

Rttlfm, p « 0.35 fm. (2.5) 



Using these basic quantities one gets from eqs. (|2.3| , |2.4|) : 



M(0) « 350 Mel/, -(^) « (250 Mel/) 3 . (2.6) 

Another quantity closely associated with chiral symmetry breaking is the pion decay constant 
which in the instanton vacuum is given by O] 



F< « 4iV c / / ^ — - r — = const" • — p z In — ~ (100 MeV) . (2.7) 
7 (2tt) 4 [M 2 (p) +p 2 ] 2 \/ F \P J 

All these quantities appear to be close to their phenomenological values. I would say 
that I don't know of any other approach to non-perturbative QCD (except, of course, brute- 
force lattice calculations) which would relate observables directly to Aqcd, arid with such an 
accuracy. Personally, I conclude that the idea that the (Euclidean) QCD partition function 
is saturated by relatively dilute instantons with quantum fluctuations of gluon field about 
them, works quite satisfactory. 
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2.2 Instanton- induced interactions 



As I have mentioned, these results have been obtained from considering the motion of light 
quarks in a given instanton background and then averaging it over the instanton ensemble 
[O. There is a mathematically equivalent technique to rederive these results, namely one 



first averages over the instanton ensemble ^M. This averaging induces many-quark interac- 
tions whose simplified version was first suggested by 't Hooft ||. Using the effective quark 
interaction theory one can calculate various observables. According to the derivation of ref. 
[0 (recently reviewed in refs. [@, 0) averaging over instanton ensemble leads to a specific 
form of the QCD partition function valid at low momenta, p < 1/p. In what follows we 
shall use the Euclidean formulation of the theory; Nf is the number of light fermion flavours 
whose masses are put to zero for simplicity. The effects of non-zero current quark masses 
have been considered in ref. |13] . 



It is convenient to decompose the 4-component Dirac bi-spinors describing quark fields 
into left- and right-handed Weyl spinors which we denote as 

^L?A)> ^L{R)faV ( 2 - 8 ) 

where / = l...Nf are flavour, a = l...N c are colour and % = 1,2 are spinor indices. Let us 
introduce the 't Hooft-like 2iVj-fermion vertices generated by J's and J's , which we denote 
by y$> , respectively. These vertices are obtained by explicit averaging over (anti) instanton 
orientation matrices Uf and over the instanton size distribution v(p). Averaging over instan- 
ton positions in d = 4 Euclidean space-time produces the overall conservation of momenta 
of quarks entering the vertex Y, hence it is convenient to write down the quark interaction 
vertex in the momentum space. There are formfactor functions F(kp) ( |2.2j ) associated with 
the Fourier transform of the fermion zero modes of one instanton, attached to each quark 
line entering the vertex. The 2iV/-fermion vertex induced by an instanton is, in momentum 
space, 

d 4 k f 4 r dH t 



r+ =Jdp u(p) JdUnlJ—^ 2npF(k f p) J -J- 2npF(l f p) 



(27T) 



■(2n) i 5(k 1 + ... + k Nf -h- ... - l Nf ) ■ U^U^e^'fe^, [i^ fafi} {k f )^ L Mf {l f )\ } ; (2.9) 

for the Y~ vertices induced by J's one has to replace left-handed Weyl spinors ?Pl,iPl by 
right-handed ones, ipn, ip^. Using these vertices one can write down the partition function to 
which QCD is reduced at low momenta, as a functional integral over quark fields [[14], |2], [13 



Z = D^Dtf exp / d 4 x V i, f i^ f — ^- — ^- (2.10) 




where N± are the number of J's and J's in the whole d = 4 volume V. The volume factors 
in the denominators arise because of averaging over individual instanton positions, and 
certain mass factors M^ f are put in to make eq. ( ^.10| ) dimensionless. Actually, the mass 
parameter Mi plays the role of separating high-frequency part of the fermion determinant 
in the instanton background from the low-frequency part considered here. Its concrete value 
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is irrelevant for the derivation of the low-energy effective action performed below; in fact 
it is established from smooth matching of high- and low-frequency contributions to the full 



fermion determinant in the instanton vacuum |T0 . 

Having fermion interactions in the pre-exponent of the partition function is not conve- 
nient: one should rather have the interactions in the exponent, together with the kinetic 
energy term. This can be achieved by rewriting eq. ( |2.1U| ) with the help of additional inte- 
gration over 'Lagrange multipliers' A±: 



Z = Jdi/jD^ exp jiV + |hi 



N, .\ „ (, N. 



AM In =^-1 + 



\+VMx* ) V \-VMi 1 



+ [d 4 xJ2 $fi^ f + KY^ + \-Y Nf . (2.11) 
/=i J 

Since N± ~ V — > oo integration over A± can be performed by the saddle-point method; the 
result is eq. ( |2.10|) we started from. 



As seen from eq. ( ^.llj ) A± plays the role of the coupling constant in the many-quark in- 
teractions. It is very important that their strength is not pre-given but is, rather, determined 
self-consistently from the fermion dynamics itself; in particlular, the saddle-point values of 
\± depend on the phase quarks assume in the instanton vacuum. As shown below, in the 
chiral symmetry broken phase the values of A±, as determined by a saddle-point equation, 
appear to be real. 

To get the 2A^j-fermion vertices (|2.9| ) in a closed form one has to explicitly integrate over 
instanton orientations in colour space. For the A^j-fermion vertex one has to average over 
Nf pairs of (U, W). In particular, one has: 



dU=l, I du u?ul j = ±s$sl, J du uzugugug 



etc. (2.12) 



iVf- 

We present below the resulting vertices for Nf = 1, 2, 3 and for any Nf but N c — > oo. 
N f = 1 



In this case the "vertex" (|2.9| ) is just a mass term for quarks 



(2.13) 



Y ^ = W c Jj^J d P "(p)i^pF(kp)} 2 [Pl(k)^r(k) , 
with a momentum dependent dynamically-generated mass M(k) given by 

M(k) = Y c J d P v ^ [^(*P)1 2 « Y e ^ 27Tf)F ^ 2 ■ ( 2 ' 14 ) 

In all our previous work on the instanton vacuum we have assumed that the distribution 
in the sizes of instantons, v(p), is a sharp function peaked at certain p, and replaced p 
by this p in the argument of the formfactor functions F(kp). However, there is a subtlety 
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here: if the size distribution for large p behaves as u(p) ~ 1/p 3 (corresponding to the linear 
potential between heavy quarks fL6| ) the dynamical quark mass logarithmically diverges at 
small momenta implying, in a sense, the confinement of light quarks. We shall not pursue 
this interesting topic here but replace each time p by its average value p. 

In order to find the overall scale A of the dynamical mass one has to put ( p. 13 ) into 



eq. fl2.11p , integrate over fermions, and find the minimum of the free energy in respect to 



A±. At 9 = the QCD vacuum is CP invariant so that N + = N_ = N/2 and consequently 
A + = A_ = A Q. In this case the 75 term in gets cancelled, and the exponent of the 
partition function ( 2.11]) reads: 



-AHnA + J 'd'x /-J^Tr In ^ + l ^pF{kp)f 
-^lnA + 2Ar c v|^ln|A; 2 +^ [ 2 vr pF(A;p)] 2 ) |. (2.15) 



Differentiating it in respect to A and using eq. (|2.14j) one gets the gap eq. ( |2.3|) or the self- 
consistency condition which is in fact a requirement on the overall scale of the constituent 
quark mass M(k); its momentum dependence is anyhow given by eq. Q2.14 ). Since the 



momentum integration in eq. fl2.3|) is well convergent and is actually cut at momenta k ~ 1 /p, 
the saddle-point value of the 'Lagrange multiplier' A is of the order of \J N C N/V /p. The 
steepness of the saddle-point integration is proportional to the volume V, hence the use of 
the saddle-point method is absolutely justified. 

Note that eq. ( |2.13p reproduces the massive quark propagator ( J2.1|) , hence the chiral 
condensate is given by eq. ( |2.4j) . It is very important that, initially, one does not know 
the strength of the quark interactions (represented by the 'Lagrange multiplier' A): it is 
fixed only after integration over the quark fields is performed. We also stress that the basic 
quantities associated with spontaneous chiral symmetry breaking, such as < ipip >,M or 
F n are non-analytic in the instanton density, N/V: such a behaviour is characteristic of 
spontaneous breaking of continuous symmetry. 

N f = 2 



In this case averaging eq. ( |2.13D over the instanton orientations with the help of eq. Q2.12D 



gives a nontrivial 4-fermion interaction. It is, of course, non-local: a formfactor function 
F(kp) is attributed to each fermion entering the vertex; in addition it should be averaged 
over the sizes of instantons. The non-locality is thus of the order of the average instanton 



size in the vacuum. One has [14, 15 



y+_ t l ^k^k 2 d%d% 



dp u(p) (2np) 4 F{k lP )F{k 2 p)F{hp)F{l 2 p) 



fluctuations of the topological charge, N + — leading to the so-called topological susceptibility 
(related to the solution of the U(l) problem) has been considered in refs.||l4|, Q. 
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2N 



8N 



(2.16) 



For the /-induced vertex Y~ one has to replace left-handed components by right-handed 
ones. In all square brackets summation over colour is understood. Note that the last-line 
(tensor) term is suppressed at large N c ; it, however, is crucial at N c = 2 to support the actual 



SU(4) chiral symmetry in that case 15|| . The antisymmetric e^^ 2 e 91S2 structure demonstrates 



that the interactions have a determinant form in the two flavours. Using the identity 



2e hh 



5 fl 5 h 

91 92 



(r A ) h (r A ) h 



(2.17) 

and adding the /-induced vertex one can rewrite the leading- N c (first) term of eq. Q2.16p 



as 



(2.18) 



which resembles closely the Vaks-Larkin [17| / Nambu-Jona-Lasinio |1| model. It should 
be stressed though that in contrast to that at hoc model the interaction ( |2.16| ) i) violates 
explicitly the ?7a(1) symmetry, ii) has a fixed interaction strength related to the density of 
instantons (see below) and Hi) contains an intrinsic ultraviolet cutoff due to the formfactor 
functions F(kp). In addition, at N c = 2 it correctly preserves the actual S*?7(4) chiral 
symmetry. We shall show in the next section that the four-fermion interaction (|2.16|) leads 
to the spontaneous chiral symmetry breaking, with the appearance of the constituent quark 
mass M(k) satisfying the same gap equation (2^) as we obtained above in the case Nf = 1. 



N f 



In this case one gets a 6-fermion vertex of the following structure [ |l5|j : 

„-3 



Y, 



d i kid i k2d i k^dHidH2d 4L U 
' = N C {N* -I) J (2^p 



5(ki + k 2 + k 3 -li-l 2 - h 



-3! e 



/1/2/3 



e 9l9293 1 1 



dp v{p) (27rp) e F(k 1 p)F(k 2 p)F(k 3 p)F(hp)F(l 2 p)F(l 3 p) 
3 



2(iV c + 2) 



[< 1 (A:iX(/i)]^i /2 (A; 2 )^ 2 (/ 2 )]^i / 3(fc3X(/3)] 



+ 



^x^ttlnihMihmU^ (2.19) 

Again, it is of a determinant structure, this time in 3 flavours, and again the tensor term is 
suppressed at N c — ► 00. 

Any N f 

For arbitrary Nf the leading term at N c —>■ 00 can be written as a determinant of Nf x Nf 



matrices composed of non-local chiral quark bilinears Jf g |L4j , |15| , [131 : 



Figure 1: Instanton- induced 2iVy-quark vertex. The black blobs denote the formfactor func- 
tions F(kp) attached to each quark leg. 



where 



Y, 



N r 



Mr 



d 4 x dpu(p) detAr iJ ± (x ) p) 



(2.20) 



d A kdH 



[27TpF(kp)}[27TpF(lp)} 



If one neglects formfactors, that is puts F = 1, these chiral currents become local, 



J±(*,p)^(27Tp) 2 



(2.21) 



(2.22) 



2.3 Bosonization 

One can linearize the many-fermion vertices induced by instantons by introducing auxiliary 
boson fields. This formal procedure is called bosonization of the theory. Roughly speaking, 
when one has a theory with 4-fermion interactions, it can be viewed as a limit of a one-boson 
exchange when the mass of intermediate boson tends to infinity. This is the meaning of the 
bosonization. 

In case Nf = 2 the instanton-induced interactions are 4-fermion ones (see eqs. fl2.16 



2.181) ) and it is very easy to bosonize them by introducing scalar and pseudoscalar fields. 
Note that the non-leading second term in eq. ( |2.16| ) requires additional tensor fields for the 
bosonization. In case Nf > 3 the instanton-induced interactions become 6-fermion and so 
forth, so that the bosonization becomes less trivial. However, it is still possible to perform 
it using N c as an algebraically large parameter [14], [15], Indeed, introducing Nf x Nf 
matrices M. the following equation becomes true in the saddle-point approximation (justified 
at large iV c ): 



exp Adet 



N. 



dM exp <^ iTi(MJ) - (N f - 1) 



'det[.MiV c 



(2.23) 
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This remarkable formula enables one to bosonize many-fermion interactions of the deter- 
minant type, like in eq. (|2.20|) . It should be stressed however that the procedure is justified 
only at large N c , otherwise i) the fermion interactions have not the simple determinant form 
and ii) the saddle-point evaluation of the integral in eq. ( 2.23|) is not justified. We notice 
that at Nf = 2 eq. ( [2.231 ) becomes exact since in this case the integral over M is Gaussian. 
Indeed, in this particular case the power of det M is unity while the determinant of a 2 x 2 
matrix itself is quadratic in matrix entries. 

We are now prepared to write down the partition function ( |2.11| ) by introducing auxiliary 
integration over 'meson' fields Ml,r(x) coupled linearly to the quark chiral currents J^x); 
both quantities are Nf x Nf matrices in flavour. We have [15, IHfl : 



Z= [— exp(-iVlnA) 
J 2tt 



jDM L , R exp Jd 4 xl-(N f -l] 



'det[A4 L iV c ] 



'det[M R N c 



D^Dip exp J d A x {^ f i$ip f + iTr [M L J + ] + iTr [M R J-]}. (2.24) 

The last line presents a theory of quarks interacting with external chiral 'meson' fields, 
M.L.R- Integration over quarks can be performed by expanding the resulting functional 
fermion determinant in powers of Ml,r and / or their derivatives. A concrete example of 
such expansion will be given below. 

The second line in eq. (|2.24|) presents the potential energy of the Ml,r fields: it has a 
rather peculiar form of a power of the determinant composed of these fields. In the particular 
case Nf = 1 these terms vanish, as it should be since at Nf = 1 instantons induce just a 
mass term for fermions, see the previous subsection. In the particular case of Nf = 2 the 
potential energy of the Ml,r fields becomes quadratic in the fields, that is to say, we have 
just a mass term for the 'meson' fields. Notice that there is no kinetic energy term at the 
tree level: the kinetic energy (as well as higher derivative terms) is generated dynamically 
after one integrates over quarks, that is through quark loops. 

The fermion action (last line in eq. ( |2.24| )) is invariant under full chiral rotations with 
arbitrary U(Nf) matrices A, B: 




4 



4^, 



Mr 



BM l A\ 
AM R Bl 



(2.25) 



However, the potential energy of the 'meson' fields (the second line in eq. ( p. 24] )) has a 
smaller invariance. Indeed the determinants transform as 



det[M L ] -> detL4 f 5] det[M L ], det[M R ] -> det[#U] det[M R ], (2.26) 

therefore they acquire a U(l) phase factor of the relative A^B transformation. For that reason 
eq. (|2.24j ) breakes explicitly the axial Ua{X) symmetry, as expected on general grounds from 
instantons ||. 

Let us show that eq. (|2.24j ) leads to the spontaneous breaking of the chiral symmetry. To 
that end, let us parametrize the 'meson' fields by 
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M L (x) = [<t(x) + T](x)]U(x)V(x), , , 

M R {x) = [a(x) -n{x)]V(x)W{x) 1 > 

where cr(x) is the scalar flavour-singlet field and 7](x) is the pseudoscalar flavour- singlet field. 
The SU(Nf) matrix fields U(x) and V(x) can be further on parametrized by scalar a (x) 
and pseudoscalar it (x) fields belonging to the adjoint representation of the flavour SU(Nf) 
group: 



U(x) = exp iir A (x)\ A 
V(x) = exp ia A (x)\ 



A f ^A (2-28) 



where X A are the Hermitean generators of the SU(Nf). 

The vacuum state of the theory given by the partition function ( p. 24 ) corresponds to a 



non-zero value of the flavour-singlet a field. This is equivalent to the spontaneous breakdown 
of chiral symmetry, and the n A fields become Goldstone particles. 

To reveal the nonzero value of the a field let us calculate the effective potential for it. 
Putting 

a = const, n = 0. U = V = 1 (2.29) 

and integrating over quarks in the constant background scalar field we obtain the effective 
potential V e ff(a, A) to which we add the 'Lagrange multiplier' piece, iVm A (see eq. ( |2.24|) ): 

V eff (a,X) = y In A + 2{N f -l){aN e )*r* -2N f N c J-0~ 4 \n{k 2 + a 2 [2ixpF{kp)} A }. 

(2.30) 

We have now to minimize (|2.30| ) in both quantities, a and A. The extremum condition gives 



the following equations for the saddle-point values ctq, Aq: 



N r d 4 k M 2 (k) 

V- 4N 'I(^WTWWy ^> = M2TO)] 2 , (2.31) 
A "f^f. (2 . 32 ) 



2V V N 

These equations demonstrate that the scalar field develops a nonzero v.e.v. (7 > quarks 
get a dynamical mass M(k) and chiral symmetry is broken. The first equation is the familiar 
self-consistency or gap equation (|2.3|) which we have also obtained in the Nf = 1 case: it 
determines the overall scale of the dynamically generated momentum-dependent constituent 
quark mass M{k). The second eq. ( |2.32| ) relates the value of the 'Lagrange multiplier' A 
to the v.e.v. of the a field. We find that these quantities have the following parametric 
dependence on the basic characteristics of the instanton vacuum, their density, N/V, and 
their average size, p: 



M(0) «^ ~,/^, Ao~^(^|V. P.33) 



VN F V 
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It is important that the strength A of the 2iV/-fermion instanton- induced interactions is 
not fixed beforehand but is, rather, determined by the phase the fermion system assumes. 
In a given phase, like the chirality broken phase, the coupling constant A is defined unambi- 
giously through the extremum conditions ( |2.31| , [2~32j ). In another phase, say, chiral invariant 
or in a phase where diquarks condense 
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the saddle-point value of A would be, gener- 
ally speaking, different. It means that the formalism presented here is different from the 
'quenched approximation': it incorporates the back reaction of fermions on the instanton 
medium. The main assumption in the starting formula for this derivation, eq. fl2.i0| ), is that 
one can average independently over the collective coordinates of J's and J's . It implies that 
correlations between pseudoparticles coming from the gauge sector are neglected or, better 
to say, treated a la variational principle resulting in an effective size distribution z/(p) || |13 



however correlations arising from fermions are taken into account. 

In getting the equation for the v.e.v. of the scalar field a we have used the mean-field 
approximation, eq. ([2.31|) . Theoretically speaking, its accuracy is of the order of 0(Nf/N c ) 
since the meson loops have been disregarded. In principle, they could be taken into account. 
It should be reminded, though, that eq. ( |2.24|) itself has been derived in the limit of N c — > oo; 
therefore, if one wishes to take into account higher order corrections in Nf/N c one should 
rather work with the unabridged vertices (|2.16|) or (|2.19|). 



2.4 Chiral lagrangian 



The low- momentum QCD partition function ( |2.24| ), after integration over quarks in the given 
background 'meson' fields M. l,r, gives an effective action for Nj scalar and Nj pseudoscalar 
fields which one can parametrize according to eqs.( 2.27 , |2.28| ). The scalar flavour-singlet field 
a(x) should be counted from its mean-field value <r given by eq. ( |2.31| ) 
multiplier' A should be put to its saddle-point value (|2.32|). 



and the 'Lagrange 



The 2Nj fields introduced by eqs.( 2.27 , |2.28| ) are not properly normalized: to get the 
correct normalization one has to extract their kinetic energy (= two derivative) terms from 
the quark functional determinant and redifine the fields so that the kinetic energy is standard, 
|(^0) 2 . The cubic, quartic,... terms in the boson fields arising from eq. ( 2.24 ) will then give 



the meson interactions. It can be shown from N c power counting []14j that the cubic coupling 
constants for the properly normalized meson fields are of the order of l/y/N C) the quartic 
couplings are ~ 1/N C , and so on. This is as it should be from the general N c counting rules. 

We have already explained in the previous subsection that the axial {/a(1) symmetry 
is explicitly (not spontaneously!) broken in eq. fl2.24|) , therefore the pseudoscalar flavour- 
singlet r) meson is not a Goldstone boson, the Ua(1) problem is solved. Moreover, in the limit 
Nf/N c — > we recover from eq. ( |2.24j ) |TjL |T3| the theoretical Witten-Veneziano formula for 
the singlet rf mass, as given by 



2N f < Q| > /V 



(2.34) 



where < Q\ > /V =< (N + — iV_) 2 > /V is the topological susceptibility of the vacuum. 



As to the non-singlet pseudoscalar fields n (x) introduced by eqs.( [2.27| , |2.28| ) which we 



shall call pions for short, they appear to be massless Goldstone fields. Indeed, the constant 



fields ti a correspond to global chiral rotations ( p.25|) with = B = exp(i7r A A A ). Eq. ( p.24|) 
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is invariant under such rotations, therefore the lagrangian for the ir A (x) fields contain only 
derivatives, i.e. pions are massless - in accordance with the Goldstone theorem. 

One can check from eq. Q2.24Q fll4|l that scalar meson fields get the mass of the order of 



1/p; numerically it is in the 1 GeV range. Though the rj mass, algebraically, is given by a 
different formula ( j2.34| ), numerically it also turns out to be about 1 GeV. Strictly speaking, 



in that range of momenta the low-energy QCD partition function ( |2.24| ) is not justifiable 
as it has been derived above from the partition function (|2.10|) valid for momenta k < 1/p. 
Therefore, in a consistent approach to the low-momentum theory at k < 1/p one has to 
freeze out all the meson fields except massless pions. [In the academic limit of Nf/N c — > 
one would also need to keep the r( degree of freedom.] It is very important that the quark 
masses are, parametrically speaking, much less than 1/p: the dimensionless quantity 

(Mp) 2 ~ ^fp'y « 1 (2.35) 

is suppressed by the packing fraction of instantons in the vacuum. The whole approach to the 
instanton vacuum implies that instantons are on the average relatively dilute and that this 
packing fraction is numerically small. Theoretically, the smallness of the parameter ( p.35|) 



can be traced back to the "accidentally" large coefficient in the Gell-Mann-Low function, 



the famous 11/3 ||, [13]. 



We, thus, arrive to the conclusion that at low momenta k < 1/p there are exactly two 
degrees of freedom left: quarks with a dynamical mass M <C 1/p and the massless Goldstone 
pions. In principle, one could think of a 'soft' instanton size distribution going at large sizes 
as z/(p) ~ 1/p 3 - Such a distribution is peculiar because it automatically leads to a linear 
potential between heavy quarks |16| and could therefore reproduce confinement. Moreover, 



there are certain reasons to believe that this particular distribution is, effectively, realized 
in nature [Diakonov and Petrov, in preparation]. As explained above in such a case M(k) 
would logarithmically diverge at k — > (see eq. (p. 14 )) so that free quarks would not show up. 



Amusingly, such a possibility would not invalidate the use of the effective chiral theory for 
bound state problems, like inside the nucleons, since for virtualities k ~ (nucleon radius) _1 
the estimate (|2.35|) would still hold true. In these lectures, however, we shall not pursue 
this interesting possibility of marrying confinement with the chiral theory but assume that 
averaging over instanton sizes merely replaces p by its peak or average value p. 

Having made these preliminary remarks, let us write down the effective partition function 
to which QCD is reduced at low momenta, k < 1/p. It follows from the instanton-induced 
partition function ( [2.24| ) where we freeze out all meson fields except pions and put a and A 
to their saddle-point values: 

Z = Jdtt a jD^D^exp Jd 4 x j ^ f (x)i^ f (x) + % J d j^e i ^ x) ^ M(k)M{l) 
4^)(^^)^ + C/ ] / ^) 1 ^)^ a (0]}, Uf(x) = (e Wi n A (x)X A )) f g . (2.36) 



This effective theory has been first derived in refs. [TT|, [14]]. Eq. ( |2.36| ) shows quarks in- 
teracting with chiral fields U(x), with formfactor functions equal to the square root of the 
dynamical quark mass attributed to each vertex where U(x) applies. The matrix entering in 
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the parentheses is actually a Nf x Nf matrix in flavour and a 4 x 4 matrix in Dirac indices. 
It can be identically rewritten as 



U(x)^^ + ?7t( x )l_T5 = exp (i-k a {x)\ a ^) = U^'(x), (2.37) 

the industrious final abbreviation being due to Pavel Pobylitsa. 

The formfactor functions y/M(kp) for each quark line attached to the chiral vertex auto- 
matically cut off momenta at k > 1/p. In the range of quark momenta k ^ 1/p (which we 
shall be mostly interested in) one can neglect this non-locality, and the partition function 
( |2.36|) is simplified to a local field theory: 



Z = J Dtt a J Dip* Dip exp J d 4 x ip\x) [i$ + iMW 5 (x)) ip{x). (2.38) 

One should remember, however, to cut the quark loop integrals at k ~ 1/p ~ 600 MeV. 
Notice that there is no kinetic energy term for pions: it appears only after one integrates over 
the quark loop, see below. Summation over colour is assumed in the exponent of eq. (|2.38f) . 

Eq. ( |2.38| ) defines a simple and elegant local field theory though it is still a highly non- 
trivial one. Its main properties will be established in the next section Q. 



3 Properties of the effective chiral lagrangian (EChL) 

Properties of effective theories of quarks interacting with various meson fields have been 
studied by several authors in the 80's, most notably by Volkov and Ebert and Dhar, 
Shankar and Wadia [ 2"2f . The fact that integrating over quarks one gets, in particular, the 
so-called Wess-Zumino term has been first established by Eides and myself [p0| |, though in 
somewhat different settings, see also below. 

Integrating over the quark fields in eq. ( 2.38 ) one gets the effective chiral lagrangian 
(EChL): 



S, 



eff 



7T 



-iV c lndet UB + iMU^ 



(3.1) 



5 We have been asked about the relation of this low-energy theory with that suggested by Manohar and 
Georgi Jig] . One can redefine the quark fields 



tp — > ip' = cxp(iir A \ A j 5 /2)i/j 
and rewrite the lagrangian in ( |2.38| ) as 

£ = ^'(ip+Y 

with 

i 

A 



4 $P = ^ exp(i7r A A A 7s/2), 
y 5 + iM)ijj' 



(2.39) 



(2.40) 



(2.41) 



which resembles closely the effective lagrangian of Manohar and Georgi (the effective chiral lagrangian in a 
similar form has been independently suggested in ref. ]20|]). 

The crucial difference is that Ma nohar and Georgi have added an explicit kinetic energy term 
i^Tr (df.Wdf.U)/^ on top of eq. (^40|) . This is a typical double counting as the kinetic energy term arises 
from quark loops, see the next section. 
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The Dirac operator entering eq. ( |3.1| ) is not hermitean: 



D = i$ + iMW 5 , D ] = i$- iMlT*\ (3.2) 
therefore the effective action has an imaginary part. The real part can be defined as 



ReS, 



eff 



k] — In det ^ 

- J 2 



D*D = -d 2 + M 2 — M(^7 75 ) 



dId = 



-d 2 + M 2 



(3.3) 



In the next two subsections we establish the properties of the real and imaginary parts 
of the EChL separately, following ref. |23]. 



3.1 Derivative expansion and interpolation formula 



There is no general expression for the functional (|3.1| ) for arbitrary pion fields. For certain 
pion fields the functional determinant ( |3.1|) can be estimated numerically, see section 4. 
However, one can make a systematic expansion of the EChL in increasing powers of the 
derivatives of the pion field, dU. It is called long wave-length or derivative expansion. 
Moreover, one can do even better and expand the real part of the EChL in powers of 



pM 



+ M 2 



(17-1) 



(3.4) 



where p is the characteristic momentum of the pion field. This quantity becomes small in 
three limiting cases: i) small pion fields, 7r (x) <C 1, with arbitrary momenta, it) arbitrary 
pion fields but with small gradients or momenta, p <C M, Hi) arbitrary pion fields and 
large momenta, p 3> M. We see thus that expanding the EChL in this parameter one gets 
accurate results in three corners of the Hilbert space of pion fields. For that reason we call 
it interpolation formula . Our experience is that its numerical accuracy is quite good for 
more or less arbitrary pion fields, even if one uses only the first term of the expansion in 



(|3.4|) , see below. 

The starting point for both expansions is the following formal manipulation with the real 
part of the EChL ( |3.3| ) . The first move is to use the well-known formula, In det [operator] = 
Sp In [operator], where Sp denotes a functional trace. One can write: 



N c 

HeS e ff[ir] = — ^ In det 



d 2 + m 2 )- 1 m(^c/ 7s ; 



-ySpln [l - (-d 2 + M 2 )- l M($U^\ 



2 J J (2vr) 4 
d^x I — , , Tr In 



-ik-x Tl ln 



-d 2 + M l y Y M{$U 



N c f j4 f #k 



1 — (Ar + M — 2ik ■ d — d^M^U 1 " 



Ak-x 



1. 



(3.5) 



2 J ~ ~ J (2tt) 4 

In going from the second to the third line we have written down explicitly what does the 
functional trace Sp mean: take matrix elements of the operator involved in a complete basis 
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(here: plane waves, exp(ik-x)), sum over all states (here: integrate over c? 4 A;/(27r) 4 ) and take 
the trace in x. 'Tr ' stands for taking not a functional but a usual matrix trace, in our case 
both in flavour and Dirac bispinor indices. In going from the third to the last line we have 
dragged the factor exp(ik ■ x) through the operator, thus shifting all differential operators 
d — > d + ik. We have put a unity at the end of the equation to stress that the operator is 
acting on unity, in particular, it does not differentiate it. The above is a standard procedure 



for dealing with functional determinants [24 



The last expression in eq. Q3.5| ) can be now easily expanded in powers of the derivatives 



of the pion field: it arises from expanding (|3.5| ) in powers of (jiU" 15 and of 2ik • d + d . The 
first non-zero term has two derivatives, 



ReS^ M - fd*x Tr ( M ^ 

„ + „ „n .„ r d A k M 2 



2 



4 J*** W>' u ) ' 4JV = S^f WTWf (3 ' 6) 

It is the kinetic energy term for the pion field or, better to say, the Weinberg chi- 
ral lagrangian: atually it contains all powers of the pion field if one substitutes U(x) = 
exp(i7r (x)A ). The proportionality coefficient (the last factor in eq. (|3.6|) ) is called F%, ex- 
perimentally, F n ~ 94MeV. The last factor in eq. ( |3.6|) is logarithmically divergent; to make 
it meaningful we have to recall that we have actually simplified the theory when writing it in 
the local form (|2.38|) . Actually, the dynamical quark mass M is momentum-dependent (see 



eq. ( p. 36 )); it cuts the logarithimically divergent integral at k ~ 1/p. Using the numerical 



values of p w 600 MeV and M w 350 MeV we find 

F 2 = AN / — — — « — M 2 ln— w (100 MeV) 2 (3 7) 

being not in a bad approximation to the experimental value of F n . Actually, the two- 
derivative term is the only divergent quantity in the EChL: higher derivative terms are all 
finite. 

A more standard way to present the two-derivative term is by using hermitean Nf x Nf 
matrices = iWd^U. One can rewrite eq. (|3.6Q as 



[tt] = S I d *xTr L^, = lU^d^U. (3. 



The next, four-derivative term in the expansion of ReS e ff is (note that the metric is 
Euclidean) 

N t 

ReSi%[7r] = Jd*x[2Tr(d lt Lj 2 + r FrL lt L v L„L v ]. (3.9) 

These terms describe, in particular, the d-wave rnr scattering lengths, and other ob- 
servables. They can be compared with the appropriate phenomenological terms in the 
Gasser-Leutwyler chiral perturbation theory [2S|: eq. ( |3.9| ) with the concrete coefficients 
like iV c /i927r 2 appears to be in good agreement with the phenomenological analysis p6 
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The next-to-next-to-leading six derivative terms following from eq. ( |3.5| ) have been com- 
puted in refs. [p7 , 28| ; however, a detailed comparison with phenomenology is still lacking 
here. 

I would like to mention an interesting paper |^9| where the derivative expansion of the 
EChL has been obtained from dual resonance models. For reasons not fully appreciated the 
dual resonance model for the nir scattering gives (numerically) very similar coefficients as 
those following from eq. (|3.9| ), however there is a discrepancy at the 6-derivative level. 

I now turn to the interpolation formula promised in the beginning of this subsection. 
One can start from the last line in eq. (|3.5| ) and expand in powers of M{<j)U lb ). It is clear 
that the actual expansion parameter will be ( |3.4|) . In the first non-zero order we get [^] 



ReS«r"°'[7r] 



— Id'x 
4 J J (2vr) 4 



1 



(k + id) 2 + M 2 



M($U 



75' 



1 



(k + id) 2 + M 2 



M($U 



75' 



(3: 



0) 

It will be convenient now to pass to the Fourier transform of the U (x) field understood as a 
matrix, 



U(p) = / d 4 x e ip ' x [U(x) - 1] . 



(3.11) 



The partial derivatives appearing in eq. ( |3.10|) act on the exponents of the Fourier transforms 
of U, W and become corresponding momenta. As a result we get 



ReS. 



Interpol 
eff 



71 



d 4 p 

W) 4 



p 2 Tr \u\p)U(p) 



■AN r 



d 4 k 



M 2 



(2tt) 4 (jfe 



f) 2 + M 2 



{k + l) 2 + M 2 



(3.12) 

At p — > the last factor becomes F 2 (see eq. (|3.7| )), and eq. (|3.12|) in nothing but the first 
term in the derivative expansion, eq. (^.6|). However, eq. ( |3.12| ) also describes correctly the 
functional S e ff[ir] for rapidly varying pion fields (with momenta p 3> M) and for small pion 
fields of any momenta, when one can anyhow expand eq. (|3.5|) in terms of tc a (x) and hence 
in U(x) — 1. The logarithmically divergent loop integral in eq. ( |3.12| ) should be regularized, 
as in eq. (|3.7| ). 

Similarly, one can get the next term in the 'interpolation' expansion which will be quartic 
in U(p), however our experience tells us that already eq. ( |3.12| ) gives a good approximation 
to the EChL for most pion fields. 



3.2 The Wess— Zumino term and the baryon number 

We now consider the imaginary part of S e ff[ir}. The first non-zero term in the derivative 
expansion of IrnS' e //[7r] is [^U], the Wess-Zumino term [P0| . It is known |JT| that it cannot 
be written as a d = 4 integral over a local expression made of the unitary U(x) matrices, 
however the variation of the Wess-Zumino term is local. For this reason let us consider the 
variation of lmS e ff[ir] in respect to the pion matrix U(x). We have [^J 



8 ImSeff [tt] = -N c 5 Im In det D = 
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iN r . 



Sp \{D^D)- l D^5D - (DD^D^SD^ 



(3.13) 



Now one can put in explicit expressions for D,D^ from eqs.( |3.2| , |3.3| ). The aim of this 
excercise is to get c/)U in the denominators so that an expansion in this quantity similar to 
that of the previous subsection could be used. 

Using the Dirac algebra relations (in Euclidean space) 



{-,,.-„//} 2d,,,, -;, : - /( . 

{75, 7m) = °> 7s = 7s, tr (757a7/37 7 75) = 4e a/37<5 , 
one gets after expanding eq. (|3.13|) in powers of ^U" 75 the first non-zero term 



7 + 



75 = 71727374 



4, 



5lmS, 



eff 



7T 



iN c 
48vr 2 



d A xe al3lS Ti (d n U f d«UdM f d*UU^5U 



(3.14) 



(3.15) 



It can be easily checked that this expression coincides with the variation of the Wess-Zumino 



term written in the form a d — 5 integral [31 | 



ImS, 



eff 



240vr 2 



d 5 x e 



Tr (rfdjf) (tfdpU) (U^d^u) (U ] 3 5 U) (Wd £ U) 



+ 



higher derivative terms. (3.16) 

In fact the integrand in eq. ( |3.16| ) is a full derivative, however, to write it explicitly one 
would need a parametrization of the unitary matrix U. The expansion of eq. Q3.16j ) starts 
from the fifth power of tt a (x), and it is non-zero only if Nf > 3 [ SD| . It is important that, 
similar to KeS e ff, the imaginary part is also an infinite series in the derivatives. 

The EChL ( |3.1| ) or, more generally, the low-energy partition function ( |2.24| ) from where 
eq. ( [3.1D has been derived, is invariant under vector flavour-singlet transformations. There- 



fore, there should be a corresponding conserved Nother baryon current, B^. This current 
is associated with the imaginary part of S e ff only; since lmS e ff is an infinite series in the 
derivatives so is the associated Noether current B^. For the first Wess-Zumino term ( [3.16D 



the corresponding charge is |31 



B = -^5 Jd^e ijk Ti(U ] d l U)(U ] d 3 U)(U ] d k U)+ higher derivative terms. (3.17) 

The explicitly written term is the winding number of the field U(x). Let me briefly explain 
this notion. 

If ^(x) — > at spatial infinity so that C/(x) — > 1 in all directions, one can say that 
the spatial infinity is just one point. Eq. (|3.17p gives then the winding number for the 
mapping of the three-dimensional sphere S 3 (to which the flat d = 3 space is topologically 
equivalent when 00 is one point) to the parameter space of the SU(Nf) group. In case 
Nf = 2 the parameter space is also S* 3 so that the mapping is S 3 1— > S 3 . The topologically 
non-equivalent mappings U(x), i.e. those which can not be continuously deformed one to 
another, are classified by their winding number, an integer analytically given by eq. ( |3.17| ). 



In case of Nf > 2 mathematicians prove that mappings are also classified by integers given 
by the same eq. (|3.17|) . 
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There exists a prejudice that the baryon number carried by quarks in the external pion 
field coincides with the winding number of that field: generally speaking it is not so because 
of the higher derivative terms omitted in eq. ( p.17 ). Only if the pion field is spatially large 



and slowly varying so that one can neglect the higher derivative terms in eq. ( |3.17| ) one can 
say that the two coincide. Otherwise, for arbitrary pion fields, the baryon number is not 
related to the winding number: the former may be zero when the latter is unity, and vice 
versa. 

To see what is going on here, let us calculate directly the baryon number carried by 
quarks in an external time-independent pion field U(x) f23fl . The definition of the baryon 
charge operator in the Minkowski space is 

B = ±- Jd 3 *^. (3.18) 

Passing to Euclidean space (which we prefer to work with since functional integrals are more 
readily defined in Euclidean) one has to make a substitution ^ — > —iip^jo — > 74, so that 

B = [d 3 *^^. (3.19) 



The baryon charge in the path integral formulation of the theory given by eq. Q2.38|) is then 



B=(B} = J tfx^W) = ~i J rf3x Tr {AlS$ + iMU^Y^x) 
= — i / d 3 x Tr (x 4 , x|— : — — — — |x 4 ,x) 

=- s /Ar^ Tr( xi^-ix> 

J J -00 Z7T UJ + Itl 

= Sp 9(-H) = number of levels with E < 0. (3.20) 

Here 



H = iAlkd k + Af T4 £rw (3.21) 

is the Dirac hamiltoniam in the external time-independent pion field U(x) and 9 is a step 
function. 

Eq. ( |3.20| ) is divergent since it sums up the baryon charge of the whole negative-energy 
Dirac continuum. This divergence can be avoided by subtracting the baryon charge of the 
free Dirac sea, i.e. with the pion field switched out, H = 747^^ + M74: 



In performing the integration over u we have closed the oj integration contour in the upper 
semiplane. Had we closed it in the lower semiplane we would obtain —Sp[6(H) — 6(Hq)] 
which is the same result since Sp[9(H) + 9(-H) - 9(H ) - 9(-H )) = 0. 
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We have thus obtained a most natural result: the baryon charge of quarks in the external 
pion field is the number of negative-energy levels of the hamiltonian (|3.21|) (the number of 
the levels of the free hamiltonian subtracted). 

One can perform the gradient expansion for the baryon number similarly to that of the 
real part of the EChL. To that end let us write 



B = - / cf x 



duo 

30 2ix 



Tr (x\ 



H 







oo^ + H 2 



CO 



where 



(3.23) 



H 2 = -d 2 k + M 2 — M-i A (d k U l5 \ 



Hi 



-dl + M 2 . 



Calculating the matrix element in the plane-wave basis one gets 



B 



cf x / — 



duo r d 3 k 



2tt J (2tt) 



Tr7 4 



7 • (d + ik) + MU^ 



(3.24) 



(d + ik) 2 + M 2 -M-i- (0ET») 



7 • (d + ik) 



uO A 



■ 1. 



{d + ik) 2 + M 2 \ " ^ 3 ' 25 ^ 

For slowly varying fields U (x) eq. ( |3.25| ) can be expanded in powers of dll 15 and d (applied 
ultimately to IT 15 ). Because of the Tr7 5 ... the first non-zero contribution arises from ex- 
panding the denominator in eq. ( |3.25| ) to the third power of 7 ■ (dW 5 ). Integrals over 00 and 
k should be explicitly performed. After some simple algebra one gets 



B = "24-2 J^xeijkTridiU^djUdkU^U) + higher derivative terms (3.26) 

coinciding with eq. ( |3.17| ) derived from the Noether current corresponding to the Wess- 
Zumino term ( glBj ) ||, ||, ||, |5[. 

It should be stressed that the baryon number carried by quarks in the background pion 
field is equal to the topological winding number of the field only if it is a slowly varying one. 
The deep reason for it is the following |23|. Imagine we start from a pion field U(x) whose 



winding number is one but whose spatial size is tending to zero. Such a field would have no 
impact on the spectrum of the Dirac hamiltonian ( |3.21|) : it would remain the same as that 



of the free hamiltonian, namely it would have the upper (E > M) and lower (E < —M) 
Dirac continua separated by the mass gap of 2M. 

We now (adiabatically) increase the spatial size of the pion field preserving its winding 
number equal to unity. Since the winding number is dimensionless this can always be done. 
At certain critical spatial size the potential well for quarks formed by the external pion field 
is wide enough so that a bound-state level emerges from the upper continuum. With the 
increase of the width of the potential well the bound-state level goes down towards the lower 
Dirac continuum. Asymptotically, as one blows up the spatial size of the pion field (always 
remaining in the winding number equal unity sector) the bound-state level travels all the 
way through the mass gap separating the two continua and joins the lower Dirac sea - this 
is a theorem proven in ref. |23|1 . At this point one would discover that there is an extra state 
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close to the lower Dirac continuum (as compared to the free, that is no-field case). Therefore, 
one would say that the baryon number is now unity, - in correspondence to eqs.( |3.17| , |3.26| ). 

In a general case, however, the baryon number of the quark system is the number of 
eigenstates of the Dirac hamiltonian fl3.21| ) one bothers to fill in. The role of the winding 
number of the background pion field is only to guarantee that, if the spatial size of the field 
is large enough, the additional bound-state level emerging from the upper continuum is a 
deep one: asymptotically it goes all the way to the lower continuum. 



4 The nucleon 

4.1 Physical motivations 

All constituent quark models start from assuming that the nearly massless light quarks of 
the QCD lagrangian obtain a non-zero dynamical quark mass M ~ 350 — 400 MeV. This is 
due to the spontaneous breaking of chiral symmetry, its microscopic driving force being, to 
my belief, instantons, as explained in section 2. Even if one does not believe in instantons 
as the microscopic mechanism of spontaneous chiral symmetry breaking one has to admit 
that once a constituent quark mass is introduced such quarks inevitably have to interact 
with Goldstone pions. The lagrangian ip(i$ — M)ip is not invariant under axial rotation 
■0 — > exp(ia A X A 75 )■?/>, it is the chiral lagrangian ( [2.381 ), 

£ = MW 5 )ip, (4.1) 

which is, since the rotation of the quark fields can be compensated in this lagrangian by 
renaming of the pion fields. 

Expanding £7 75 = exp(in A \ A j 5 / F T ) in powers of the properly normalized pion fields (that 
is why we have inserted the F n constant) we see that the dimensionless constant of the linear 
coupling of pions to constituent quarks is, numerically, quite large: 

M 

g^ qq = — ~ 4. (4.2) 

I would like to emphasize that this is a model-independent consequence of saying that quarks 
get a constituent mass. 

Actually, the coupling (|4.2j ) is so strong that one may wonder how some people manage 
to get along without taking it into account. Not surprisingly, baryon models which do take 
into account pion exchange between constituent quarks give much more realistic predictions 
than, say, the old simple-minded Isgur-Karl model (for a review see ref . |36| ) . 



Moreover, at distances between quarks of the order of 0.5 fm typical for interquark 
distances inside nucleons, neither the one-gluon exchange nor the supposed linear potential 
are as large as the chiral forces. Therefore, it is worthwhile to investigate whether the chiral 
forces alone are able to bind the constituent quarks inside nucleons. Such approach may 
or may not be successful for describing high nucleon excitations where, according to the 
standard logic, the confining forces become crucial. To that we can remark the following: 
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Figure 2: Quarks in a nucleon interacting through chiral fields. 



1. The constituent quark mass is momentum- dependent; the behaviour of M(k) at low 
virtualities k may well be divergent (see section 2) and it may thus play the role of confining 
forces. [To my knowledge, this line of thought has not been pursued in the literature.] 

2. Highly excited baryons with large angular momenta, if understood as chiral solitons, 
lie on linear Regge trajectories with a realistic slope related to the F n constant [|l], j37fl . 

Therefore, it may be expected that even highly excited baryons can be incorporated into 
the chiral theory. After all, as emphasized by Witten |31| , the theory of all hadrons can 
be, in principle, formulated completely in terms of the EChL. For example, there should 
exist an EChL corresponding to the Lovelace-Shapiro dual resonance amplitudes for tttt 
scattering exhibiting the correct Regge behaviour ||28|| . A (unsolved) problem is to formulate 
an appropriate EChL in the field-theoretic language. 

Leaving aside these interesting problems, we concentrate on the lowest state with baryon 
number one, i.e. the nucleon. As mentioned above the interquark separations in the ground- 
state nucleon are moderate (order of 0.5 fm) and it is worthwhile asking whether the simple 
EChL ( |2.38| ) is capable of explaining the basic properties of the ground-state nucleon. Notice 
that the expected typical momenta of quarks inside the nucleon are of the order of M « 
350 — 400 MeV, that is perfectly inside the domain of applicability of the low- momentum 
effective theory (|2.38| ), according to its derivation in section 2. 



The chiral interactions of constituent quarks in the 3-quark nucleon, as induced by the 
effective theory ( |2.38|) , are schematically shown in Fig. 2, where quarks are denoted by solid 
and pions by dashed lines. Notice that, since there is no tree-level kinetic energy for pions 
in eq. (|2.38|) , the pion propagation occurs only through quark loops. Quark loops induce 
also many-quark interactions indicated in Fig. 2 as well. We see that the emerging picture is, 
unfortunately, rather far from a simple one-pion exchange between the constituent quarks: 
the non-linear effects in the pion field are not at all suppressed. 

At this point one may wonder: isn't the resulting theory as complicated as the original 
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QCD itself? The answer is no, the effective low-energy theory is an enormous simplification 
as compared to the original quark-gluon theory, because it deals with adequate degrees of 
freedom. Let us imagine that one would like to describe 'low-energy' properties of solid 
states, superconductivity for example. Would working with the underlying theory (QED) be 
helpful? Not at all. We know that the microscopic theory leads, under certain conditions, to 
the rearrangement of atoms into a lattice, so that translational symmetry is spontaneously 
broken. As a result the Goldstone bosons appear (here: phonons), and electrons get a 
dynamical mass different from the input one. The most important forces are due to phonon 
exchange between electrons: in fact they are driving superconductivity in the BCS theory. 
Playing with this analogy, nucleon is like a polaron (a bound state of electrons in the phonon 
field), rather than a positronium state in the vacuum. After chiral symmetry is broken we 
deal with a 'metal' phase rather than with the vacuum one, and one has to use adequate 
degrees of freedom to face this new situation. 

The instanton vacuum plays the role of the bridge between the microscopic theory (QCD) 
and the low-energy theory where one neglects all degrees of freedom except the Goldstone 
bosons and fermions with the dynamically-generated mass. Instantons do the most difficult 
part of the job: they explain why atoms in metals are arranged into a lattice and what is the 
effective mass of the electron and what is the strength of the electron-phonon interactions. 
However, one can take an agnostic stand and say: I don't care how 350 MeV is obtained 
from the microscopic Aqcd and why do atoms form a lattice in the metals: I just know it 
happens. To such a person I would advise to take the low-energy theory (|2.38| ) at face value 
and proceed to the nucleon. 

A considerable technical simplification is achieved in the limit of large N c . For N c colours 
the number of constituent quarks in a baryon is N c and all quark loop contributions are also 
proportional to N c , see section 3. Therefore, at large N c one can speak about a classical 
self- consistent pion field inside the nucleon: quantum fluctuations about the classical self- 
consistent field will be suppressed as 1 /N c . The problem of summing up all diagrams of the 
type shown in Fig. 2 is thus reduced to finding a classical pion field pulling N c massive quarks 
together to form a bound state. 



4.2 Nucleon mass: a functional of the pion field 

Let us imagine that there is a classical time-independent pion field which is strong and 
spatially wide enough to make a bound-state level of the Dirac hamiltonian ( |3.21|) for massive 
quarks, call its energy £7i eve i. We fill in this level by N c quarks in the antisymmetric state 
in colour, thus obtaining a baryon number one state, as compared to the vacuum. The 
interactions with the background chiral field are, naturally, colour-blind, so one can put N c 
quarks on the same level; the fact that one has to put them in an antisymmetric state in 
colour, i.e. in a colour- singlet state, follows from Fermi statistics. 

One has to pay for the creation of this trial pion field, however. Call this energy -Egeid- 
Since there are no direct terms depending on the pion field in the low-momentum theory 
the only origin of -Egeid is the fermion determinant ( |3.1| ) which should be calculated 



for time-independent field U(x). It can be worked out with a slight modification of section 
3. We have J53|: 
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S, jM • -X In d<A 



= -N c Sp [ln(id t + iH) - \n{id t + iH )] 

= —TN C /^Sp [ln(u; + iH) - \n(tu + iH )) , (4.3) 

where H is the Dirac hamiltonian ( |3.21| ) in the stationary pion field, H is the free hamilto- 
nian and T is the (infinite) time of observation. Using an important relation 

Sp(#-# ) = (4.4) 

(telling us that the sum of all energies, with their signs, of the Dirac hamiltonian ( |3.21| ) is 
the same as for the free case) one can integrate in eq. ( f4.3[ ) by parts and get 

duo , 



S eff [n] =TE &cld = TN c /-^Sp 

J 271 



to uo 



LO + iH lo + iHo 

= TN C J2 {En-E^). (4.5) 

£< 0) <o 

Going from the first to the second line we have closed the uj integration contour in the upper 
semiplane; owing to the trace relation ( |4.4| ) closing it in the lower semiplane would produce 
the same result. 

We see that the energy cost -Egeid one pays for a creation of the time-independent pion 
field coincides with the aggregate energy of the lower Dirac continuum in that field. The 
energy of the additional level emerging from the upper continuum, which one has to fill in 
to get the baryon number one state, .EWei, should be added to get the total nucleon mass. 
This simple scheme ^] is depicted in Fig.3. Naturally, the mass of the nucleon should 
be counted from the vacuum state corresponding to the filled levels of the free lower Dirac 
continuum. Therefore, the (divergent) aggregate energy of the free continuum should be 
subtracted, as in eq. ( |4.5| ). 

We have thus for the nucleon mass: 

M N = min (iV c £ level [7r] + £ fie i d [7r]) . (4.6) 

{tt a (x)} 

Both quantities, -Ei eve i and -Efieid, are functionals of the trial pion field 7r A (x). The clas- 
sical self-consistent pion field is obtained from minimizing the nucleon mass ( |4.6|) in 7r A (x). 
It is called the soliton of the non-linear functional ( |4.b] ), hence the Chiral Quark-Soliton 
Model. An accurate derivation of eq. ( |4.6[ ) from the path-integral representation for the 
nucleon-current correlation function is presented in ref. p3| |. It solves the problem of sum- 



ming up all diagrams of the type shown in Fig. 2 in the large- N c limit. 



6 It follows from taking the matrix trace of the hamiltonian (3.21). It should be kept in mind, though, 
that such a naive derivation can be potentially dangerous because of anomalies in infinite sums over levels. 
However, it can be checked that in this particular case there are no anomalies and the naive derivation is 
correct. 
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Figure 3: Spectrum of the Dirac hamiltonian in trial pion field. The solid lines present 
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The idea that a sigma model with pions coupled directly to constituent quarks, can be 



used to build the nucleon soliton has been first suggested by Kahana, Ripka and Soni [39 



and independently by Birse and Banerjee [(40 ] . I would call them the authors of the Chiral 
Quark-Soliton Model. Technically, however, in these refs. an additional ad hoc kinetic 
energy term for pion fields has been used, leading to a vacuum instability paradox. The 
present formulation of the model has been given in ref. [38], together with the discussion of 
its domain of applicability and its physical contents. A detailed theory based on the path 
integral approach paving the way for calcultaing nucleon observables has been presented in 
ref. 03|. 



4.3 Nucleon profile 

To find the classical pion field minimizing the nucleon mass ( ]4.6| ) one has first of all decide 
on the symmetry of the pion field. Had the field been a singlet one would take a spherically- 
symmetric ansatz. However, the pion field has flavour indices A = 1, Nj — 1. At Nf = 2 
the three components of the pion field can be married with the three space axes. This is called 
the hedgehog ansatz; it is the minimal generalization of spherical symmetry to incorporate 
the tt ± = (vr 1 ± in 2 )/^2 and vr° = tt 3 fields: 

vr^x) = n A P(r), n A = —, r= |x|, (4.7) 

r 

where P(r) is called the soliton profile. 

The choice of the ansatz is not at all innocent: baryons corresponding to different choices 
would have qualitatively different properties, see the next subsection. The maximally- 
symmetric ansatz (|4.7p will have definite consequences in applications. 
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In the Nf = 3 case there are 8 components of the 'pion' field, and there are several 
possibilities to marry them to the space axes. The commonly used ansatz is the 'left upper 
corner' ansatz Ell: 



C/(x) = exp (in A (x)y 



exp [i(n ■ r)P(r)] 







1 



n 



x 

r 



(4.8) 



As we shall see in the next subsection, the quantization of rotations for this ansatz leads to 
the correct spectrum of the lowest baryons. 

Another SU(3) ansatz [12"] discussed in the literature is (f,g,h = 1,2,3 are the flavour 
indices): 



U 



fa 



3 iP 2 /3 



cos Pi 8f g + (e tP2 — cos Pi J rifTig + sin Pi tf g hn h 



(4.9) 



where P\p{f) are spherically-symmetric profile functions. This ansatz is used to describe 
strangeness -2 dibaryons EJ[. We shall not consider it here but concentrate on the usual 
baryons for which the hedgehog ansatz ( |4.7| ) or (|4.8|) is appropriate. 

Let us first discuss restrictions on the best profile function P(r) which should mimimize 
the nucleon mass (|4.6|). What is the asymptotics of P(r) at large r? To answer this question 



one has to know the behaviour of PfieidK] for slowly varying pion fields. Using eq. fl4.5|) as 
a starting point one can work out the derivative expansion of the functional -EWd [vr] similar 



to that for the full EChL, see subsection 3.1. We have 
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+ 



2Ti{d i L l f + TiL i L j L l L j 
Li = iU^diU. 



higher derivative terms, L{ = ill^diU. (4-10) 

Substituting here the hedgehog ansatz one gets a functional of the profile function P(r); 
varying it one finds the Euler-Lagrange equation valid for slowly varying profiles, in partic- 
ular, for the tail of P(r) at large r. It follows from this equation that P(r) = A/r 2 at large 
r. The second contribution to the nucleon mass, -Ei eve i[7r], does not alter this derivation since 
the bound-state wave function has an exponential, not power behaviour at large r. Actually 
we get the pion tail inside the nucleon, and the constant A is related to the nucleon axial 



constant. This relation is identical to the one found in the Skyrme model [H: 



9a 



-AF 2 . 
3 



(4.11) 



The exponentially decreasing wave function of the bound-state level does not change this 
derivation, as well as the Goldberger-Treiman relation for the pion-nucleon coupling con- 
stant, 



9-kNN 



9aM 



N 



8ttAM 



N 



(4.12) 



Furthermore, it follows from the next four-derivative term in eq. ( f4.10| ) that the 1/r 4 cor- 



rection to P(r) at large r is absent |38|, |23[ ! It means probably that the pion tail inside 
nucleon is unperturbed to rather short distances. 
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The second important question is what should we choose for P(0)? As explained in 
subsection 3.2, a quantity which guarantees a deeply-bound state in the background pion 
field is the winding number of the field, eq. (|3.17|) . Substituting the hedgehog ansatz into 
eq. ( |3 . 1 7| ) one gets 



wind 



-- / dr snr P(r) — 
tt Jo dr 



TT 



P(r) 



sin 2P(r) 



(4.13) 



Jo 

Since P(oo) = the way to make this quantity unity is to choose P(0) = tt = 3.14.... 

An example of a one-parameter variational function satisfying the above requirements is 



2 arctg 



A = 1ri 



(4.14) 



The Dirac hamiltonian Q3.21| ) in the hedgehog pion field ( |4.7| ) commutes neither with the 
isospin operator T nor with the total angular momentum J = L + S but only with their sum 
K = T + J called the 'grand spin'. The eigenvalue Dirac equations for given value of K 2 , K 3 
have been derived in refs. [p8 , |23|j . Generally speaking, there appears a bound-state level 
with the K p = + quantum numbers whose energy can be found from solving the Dirac 
equations for two spherically-symmetric functions j, h 



dh 

— = -Mh sin P + (P levcl + M cos P) j, 
dr 



J- + _ j = Mj sin P 

dr r 



-E 



level 



M cos P) h 



(4.15) 
(4.16) 



with the boundary conditions h(0) = 1, j(0) = Cr, h(oo) = j(oo) = 0. 

These equations determine one of the two contributions to the nucleon mass, E\ eve \. The 
second contribution, namely that of the aggregate energy of the lower Dirac continuum in the 
trial pion field, which we have called Pfi e id ; can be found in several different ways. One way 
is to find the phase shifts in the lower continuum, arising from solving the Dirac equation 
for definite grand spin K |23, |45|. Another method is to diagonalize the Dirac hamiltonian 



( |3.21| ) in the so-called Kahana-Ripka basis [^6[ written for a finite-volume spherical box. 
Both methods are, numerically, rather involved. There exists a third (approximate) method 
[ P^ , P3[ ] allowing one to make an estimate of Pgeid in a few minutes on a PC. It is based on 
the interpolation formula for the EChL, see eq. ( |3.12j ) and ref. [S3] for details. 

Let us discuss the qualitative behaviour of Pi CV ei and Pfi id with the soliton scale parameter 
ro assuming for definiteness that the profile is given by eq. ( [1.14j ). 

The trial pion field plays the role of the (relativistic) potential well for massive quarks. 
The 'depth' of this potential well is fixed by the condition P(0) = tt and cannot be made 
infinite: this is related to the fact that the pion field has the meaning of angles. The 
spatial size of the trial pion field ro plays the role of the 'width' of the potential well. It 
is well known that in three dimensions the condition for the appearance of a bound state 
is MVr^ > const where V is the depth of well and 'const' is a numerical constant of the 
order of unity depending on the concrete shape of the potential well. In our case V ~ M, 
so the condition that the bound state appears is Mr ~ 1. Therefore, at small sizes r there 
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Figure 4: The classical nucleoli mass and its constituents as function of the soliton size. The 
short-dash line shows 3-Ei eve i, the long-dash line shows i£a e id, the solid line is their sum, A4n- 



is no bound state for the Dirac hamiltonian (|3.21 ), so that -Ei eve i coincides with the border 
of the upper Dirac continuum, -Ei eve i = +M. At certain critical value of Tq a weakly bound 
state emerges from the upper continuum. [For the concrete ansatz ( |4.14 ) the threshold value 
is r M ~ O.5.] As one increases r the bound state goes deeper and -Ei ev ei monotonously 
decreases. At very large spatial sizes, r — > oo, E\ cvc \ approaches the lower continuum, its 
difference from — M falling as 1/Vq [23|. The behaviour of £a eV ei as function of tq is plotted 
in Fig.4. 

The monotonous decrease of E\ evel with the increase of r is a prerogative of the trial 
pion field with winding number 1. Had it been zero, -Ei eve i would first go down and then 
start to go up, asymptotically joining back the upper continuum. In the case of N win d = — 1 
the bound state would travel in the opposite direction: from the lower towards the upper 
continuum. At N W i n d = n as much as n levels would emerge, one by one, from the upper 
continuum and travel all the way through the mass gap towards the lower one. For the trial 
pion field of the hedgehog form all these things happen exclusively for states with grand spin 

Turning now to -Efieid we first notice that for large spatial sizes of the trial pion field one 
can use the first term in the derivative expansion for -Efi e id, see eq. ( |4.10|) . On dimension 
grounds one immediatelly concludes that -EWd ~ F^ r o f° r large r , i.e. is infinitely linearly 



rising in r . At small r a slightly more complicated analysis |23 shows that -Efieid ~ r o- ^ n 



the whole, -Efi e id is a monotonously rising function of ro shown in Fig.4. 

The nucleon mass, M.n = 3E\ eve \ + -Egeid (for N c = 3) is also plotted in Fig.4 taken from 
refs. [[38, 23]. One observes a non-trivial minimum for M.^ corresponding to ro — 0.98/M ~ 
0.57 fm. This is, phenomenologically, a very reasonable value, since from eqs.( [4.11] , |4.12[ ) 
one immediatelly gets qa — 1-15 versus 1.25 (exp.) and g n NN — 13.6 versus 13.5 (exp.). The 
nucleon mass appears to be M N ~ 1100 MeV with £i evel ~ 123 MeV, E &eld ~ 730 MeV. 
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Note that the 'valence' quarks (sitting on the bound-state level) come out to be very strongly 
bound: their wave function falls off as exp(— r/0.6 fm), and about 2/3 of the quark mass 
M ~ 350 MeV is eaten up by interactions with the classical pion field. Relativistic effects 
are thus essential. 

Though the nucleon bound state appears to be somewhat higher than the free-quark 
threshold, 3M ~ 1050 MeV, there are several known corrections to it which all seem to be 
negative. The largest correction to the nucleon mass is due to taking into account explicitly 
the one-gluon exchange between both 'valence' and 'sea' quarks; this correction is 0(N C ) as 
is the nucleon mass itself. Numerically, it turns out to be about —200 MeV [|7j and seems 
to move the nucleon mass just into the right place []. 

We see thus that the 'valence' quarks in the nucleon get bound by a self-consistent pion 
field whose energy is given just by the aggregate energy of the negative Dirac continuum 
distorted by the presence of the external field. This picture of the nucleon interpolates 
between the old non-relativistic quark models (which would correspond to a shallow bound- 
state level and an undistorted negative continuum) and the Skyrme model (which would 
correspond to a spatially very large pion soliton so that the bound-state level would get 
close to the lower continuum and the field energy -Egeid would be given just by a couple of 
terms in its derivative expansion). The reality is somewhere in between: the bound-state 
level is a deep one but not as deep as to say that all the physics is in the lower Dirac 
continuum. 

Idelogically, this picture of the nucleon at large N c is somewhat similar to the Thomas- 
Fermi picture of the atom at large Z. In that case quantum fluctuations of the self-consistent 
electrostatic field binding the electrons are suppressed by large Z, however corrections go as 
powers of Z~ 2 / 3 . Therefore, the Chiral Quark-Soliton Model is in a slightly better position 
in respect to quantum corections than the Thomas-Fermi approximation. 

Apart from using the large N c approximation (which is in fact just a technical device 
needed to justify the use of the classical pion field) the Chiral Quark-Soliton Model makes 
use of the small algebraic parameter (Mp) 2 where p is the average size of instantons in the 
vacuum. This p is, roughly, the size of the constituent quark, while the size of the nucleon is, 
parametrically, 1/M. The fact that the constituent quark picture works so well in the whole 
hadron physics finds its explanation in this small numerical parameter being due to the 
relative diluteness of the instanton vacuum, which in its turn is related to the 'accidentally' 
large number (11/3) in the asymptotic freedom law [0, |I3[ . The small parameter (Mp) 2 <C 1 
makes it possible to use only quarks with dynamically generated mass and chiral fields as 
the only essential degrees of freedom in the range of momenta k ~ M <C 1/p, and that is 
exactly the range of interest in the nucleon binding problem. 



The above numerics have been obtained from the interpolation formula for |38|, ^3 



Exact calculations of -Egeid performed in |EJ as well as taking more involved profiles with 



three variational parameteres did not lead to any significant changes in the numerics. 



Following refs. |38|, |23[ there had been many calculations of the nucleon mass and of 
the 'best' profile using various regularization schemes and parameteres of the chiral model, 
|8f for a review. The effective theory derived from the instanton vacuum comes with 



sec 



7 There exist also numerous quantum corrections to the nucleon mass of different origin, which are of the 
order of 0(N®). Unfortunately, it is difficult today to treat them in a systematic fashion; see, however, the 
next subsection. 
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an intrinsic ultraviolet cutoff, in the form of a momentum dependence of the constituent 
quark mass, M(k). It can be shown on general grounds that this is a rapidly falling function 
at momenta of the order of the inverse average instanton size, 1/p ~ 600 MeV. However, 
the present 'state of the art' does not allow one to determine this function accurately at 
all values of momenta - to do so, one would need a very detailed understanding of the 
instanton vacuum. This places certain restrictions on the kinds of quantities which can 
sensibly be computed using the effective theory. Those are either finite ones, which do not 
require an UV cutoff at all, or quantities at most logarithmically divergent. Both type of 
quantities are dominated by momenta much smaller than the UV cutoff, k <C 1/p, so one 
can compute them mimicking the fall-off of M(k) by an external UV cutoff A ~ 1/p, using 
some regularization scheme. Fortunately, almost all nucleon observables belong to these two 
classes. The uncertainty related to the details of the ultra-violet regularization leads to a 
15-20% numerical uncertainty of the results, and that is the expected accuracy of the model 
today. 



4.4 Quantum numbers of baryons 

The picture of the nucleon outlined in the previous subsection is "classical": the quantum 
fluctuations of the self-consistent pion field binding iV c quarks are totally ignored. Among all 
possible quantum corrections to the nucleon mass a special role belongs to the zero modes. 
Fluctuations of the pion field in the direction of the zero modes cannot be considered small, 
and one has to treat them exactly. Zero modes are always related to continuous symmetries 
of the problem at hand. In our case there are 3 zero translational modes and a certain 
number of zero rotational modes. The latter determine the quantum numbers of baryons; it 
is here that the hedgehog (or whatever) symmetry of the ansatz taken for the self-consistent 
pion field becomes crucial. 

A general statement is that if the chiral field U c \(x) minimizes the nucleon mass functional 
( f4.6| ), a field corresponding to rotated spatial — > OijXj, or to a unitary- rotated matrix 
in flavour space, U c \ —>■ RU C \R\ has obviously the same classical mass. This is because the 
functional ( |4.6|) to be minimized is isotropic both in flavour and ordinary spaces. 

Specifically for the hedgehog ansatz [see eq. ( |4.7p for the flavour SU (2) and eq. ( f4.8|) for 



the SU(3)] any spatial rotation is equivalent to a flavour rotation. We show it for a more 
complicated case of SU(3). Indeed, the space- rotating 3x3 matrix Oy can be written as 

O^-TrO^r,) (4.17) 

where 5* is an SU(2) 2x2 matrix and Ti are the three Pauli matrices. One can immediatelly 
check that are real orthogonal 3-parameter matrices with OijOkj = Sik and OijO^ = 5jk, 
as it should be. 

When one rotates the space putting = Oij-rij the 2x2 matrix standing in the left 
upper corner of the ansatz ft4.8|) can be written as 



exp [i(n' ■ r)P(r)] = cos P(r) + i(n' ■ r) sin P(r) 

= S [cos P{r) +i(n-r) sin P(r)}Sl (4.18) 
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Therefore, if one consideres the hedgehog ansatz ( f4.8| ) rotated both in flavour and usual 
spaces, the latter can be completely absorbed into the former one: 



RU d (0x.) = RU cl (x) & 



(4.19) 



with 



R = R 



V 



s 








1 



(4.20) 



For that reason it is sufficient to consider rotations only in the flavour space. Hence there 
are 3 zero rotational modes in the SU(2) and 8-1=7 in the SU(3) flavour case. The rotation 
of the form R = exp(ia\ 8 ) commutes with the left-upper-corner ansatz and therefore does 
not correspond to any zero mode. This will have important consequences in getting the 
correct spectrum of hyperons. 

The general strategy [p3| is to consider a slowly rotating ansatz 



U(x,t)=R(t)U cl (x) R\t) 



(4.21) 



and to expand the energy of the bound-state level and of the negative Dirac continuum in 
'right' (^a) and 'left' (^a) angular velocities 

n A = -iTr (R^RX A ), tt A = -iTr (RR^X A ), fi 2 = tt 2 = 2Tr $R. (4.22) 
Taking into account only the lowest terms in the time derivatives of the rotation matrix 



R(t) one gets ]23L 49] the following form of the rotation lagrangian: 



rot 



Here Iab is the SU(3) tensor of the moments of inertia, 



Iab 



2tt 



Tr 



u + iH u + iH 



-x 1 



(4.23) 



(4.24) 



where the u integration contour should be drawn above the bound-state energy -Ei eve i to 
incorporate the 'valence' quarks. 

The appearance of a linear term in Q s is an important consequence of the presence of an 
extra bound-state level emerging from the upper Dirac continuum, which fixes the baryon 
charge to be unity. We remind the reader that in the Skyrme model this linear term arises 



from the Wess-Zumino term [|50|. For simplicity we have written unregularized moments of 
inertia, though eq. (|4.24j ) should be regularized in some way, see e.g. H49"| . 

Owing to the left-upper- corner ansatz for the static soliton (being essentially SU(2)) the 
tensor Iab is diagonal and depends on two moments of inertia, Ii 2 ' 



l AB 




A, B 
A, B - 
A, B 



-- 1,2,3, 
4,5,6,7, 



(4.25) 
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Therefore, the rotational lagrangian ( 4.23|) can be rewritten as 



rot ^Tn2 , ^ ^ o2 N c 



,E^ + yE^" 777^- ( 4 - 26 ) 

^ A=l 1 A=A 2 V3 

To quantize this rotational Lagrangian one can use the canonical quantization procedure, 



same as in the Skyrme model pO, [32l [51], |52|, E3|. Introducing eight angular momenta 



canonically conjugate to 'right' angular velocities Qa, 



and writing the hamiltonian as 

one gets 



r) T rot 



H rot = n A J A - L r (4.28) 



H rot = ± £ J\ + ^ £ 4 (4.29) 



with the additional quantization prescription following from eq. ( |4.27|) 



8 2^3 2 v ; 

In the Skyrme model this quantization rule follows from the Wess-Zumino term. In our 
approach it arises from filling in the bound-state level, i.e. from the 'valence' quarks. It is 
known to lead to the selection rule: not all possible spin and SU(3) multiplets are allowed 
as rotational excitations of the SU(2) hedgehog. Eq. ( |4.30|) means that only those SU(3) 
multiplets are allowed which contain particles with hypercharge Y — 1; if the number of 
particles with Y — 1 is denoted as 2 J + 1, the spin of the allowed SU(3) multiplet is equal 
to J. 

Therefore, the lowest allowed 577(3) multiplets are: 

• octet with spin 1/2 (since there are two baryons in the octet with Y — 1, the N) 

• decuplet with spin 3/2 (since there are four baryons in the decuplet with Y — 1, the 
A) ' 

• anti-decuplet with spin 1/2 (since there are two baryons in the anti-decuplet with 
Y = 1, the N*) 

The next are 27-plets with spin 1/2 and 3/2 but we do not consider them here. 

We see that the lowest two rotational excitations are exactly the lowest baryon multiplets 
existing in reality. The third predicted multiplet, the anti-decuplet, contains exotic baryons 
which cannot be made of three quarks, most notably an exotic Z + baryon having spin 1/2, 
isospin and strangeness +1. A detailed study of the anti-decuplet performed recently in 
j4| predicts that such a baryon can have a mass as low as 1530 MeV and be very narrow. 



Several experimental searches of this exotic baryon are now under way. 

It is easy to derive the splittings between the centers of the multiplets listed above. For 
the representation (p, q) of the SU{3) group one has 
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jlJ 2 A = \[v 2 + q 2 +vq + 3(p + ?)], (4.31) 

therefore the eigenvalues of the rotational hamiltonian (|4.29|) ) are 

= e^lf 2 + « 2 + « + 3 <f + «)] + (jJT - 2^) J < J + D - m < 4 ' 32 ' 

We have the following three lowest rotational excitations: 

(4.33) 
(4.34) 
(4.35) 

The splittings between the centers of these multiplets are determined by the moments of 
inertia, 7i )2 : 

Aio-s = £(3,o) - E&) = ^ (4-36) 

3 

A To-s = E (o% ~ E [°fi) = • (4-37) 

The appropriate rotational wave functions describing members of these multiplets are 
given by Wigner finite-rotation functions D 8,10,10 (R) H3, |5"4] . 



(p, q) = 




J = 


1/2: 


octet, spin 1/2, 


(p, q) = 


(3,0), 


J = 


3/2 : 


decuplet, spin 3/2, 


(p, q) = 


(0,3), 


J = 


1/2: 


anti-decuplet, spin 1/2 



When dealing with the flavour SU(3) case neglecting the strange quark mass m s is an 
oversimplification. In fact, it is easy to incorporate m s ^ in the first order. As a result 
one gets very reasonable splittings inside the SU (3) multiplets, as well as mass corrections 
to different observables g, g|, Q . 

In general, the idea that all light baryons are rotational excitations of one object, the 
'classical' nucleon, leads to numerous relations between properties of the members of octet 
and decuplet, which follow purely from symmetry considerations and which are all satisfied 
up to a few percent in nature. The SU(3) symmetry by itself says nothing about the relation 
between different multiplets, of course. Probably the most spectacular is the Guadagnini 



formula [50] which relates splittings inside the decuplet with those in the octet, 



8(ms* + ttin) + 3ms = 11mA + 8ms*, (4.38) 
which is satisfied with better than one-percent accuracy! 

N f = 2 case 

If one is interested in baryons predominantly 'made of u, d quarks, the flavour group is 
SU (2) and the quantization of rotations is more simple. 
In this case the rotational lagrangian is just 

^ = yE^ 2 = |E^, (4-39) 

Z i=l Z A=l 
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where the 'right' and 'left' (£Ia) angular velocities are defined by eq. ( ^4.22| ). This is 
the lagrangian for the spherical top: the two sets of angular velocities have the meaning of 
those in the 'lab frame' and 'body fixed frame'. The quantization of the spherical top is 
well known from quantum mechanics. One has to introduce two sets of angular momenta, 
Si (canonically conjugate to Qi) and Ta (conjugate to CIa)- Both sets of operators act on 
the coordinates of the spherical top, say, the Euler angles. It will be more convenient for us 
to say that the coordinates of the spherical top are just the entries of the unitary matrix R 
defining its finite-angle rotation |23[ . 

The angular momenta operators Si, Ta act on R as generators of right (left) multiplica- 
tion, 



e i(aS) j^ e -i(aS) _ e i(aa) 
e i(aS) Re -i(aS) = e -i(ar) ^ 

and satisfy the commutation relations 

[Ta,T b ] = ieABcTc, [Si, S : 

[T A ,3i] = Q, {T A , = 
A realization of these operators is 





(4.40) 
(4.41) 



2 = (^) 2 . 



(4.42) 



Si — R 



VI 



T A 



t) 

^ / pk 



kq dR 

Rkq 



pq 





dR 



(4.43) 



pq 



The rotational hamiltonian is 

H rot = QiSi - U ot 



VL a Ta — L r ' 



rot 



T 2 



(4.44) 



si _ _ 

Comparing the definition of the generators ([4.40 , 4.41 ) with the ansatz ( 4.21Q we see that 
Ta is the flavour (here: isospin) operator and Si is the spin operator, since the former acts 
to the left from R and the latter acts to the right. 

The normalized eigenfunctions of the mutually commuting operators 5*3, T3 and 5* 2 = T 2 
with eigenvalues S3, T 3 and S(S + 1) 



T(T + 1) are |3 



(S=T), 



, T353 >(R) = v^sTT(-i) T+T ^ ( SS(^) 

where D(R) are Wigner finite-rotation matrices. For example, in the S 
sentation D^ 2 (R) = R pq , i.e. coincides with the unitary matrix R itself. 
The rotational energy is thus 



T 



E 



rot 



S{S+l) T(T + 1) 



(4.45) 
1/2 repre- 



(4.46) 



2/i 2/i 

and is (2S + l) 2 = (2T + l) 2 -fold degenerate. The wave functions ( |4.45| ) describe at S — 
T = 1/2 four nucleon states (proton, neutron, spin up, spin down) and at5* = T = 3/2 the 
sixteen A-resonance states, the splitting between them being 
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m A -m N = ^- = 0(N- 1 ) (4.47) 

(coinciding in fact with the splitting between the centers of decuplet and octet in the more 
general SU(3) case, see eq. ( |4.36| )). 

It is remarkable that the nucleon and its lowest escitation, the A, fits into this spin- 
equal-isospin scheme, following from the quantization of the hedgehog rotation. Moreover, 
since N and A are, in this approach, just different rotational states of the same object, the 
'classical nucleon', there are certain relations between their properties. These relations are 
identical to those found first in the Skyrme model |44]] since they follow from symmetry 
considerations only and do not depend on concrete dynamics which is of course different 
in the naive Skyrme model. For example, one gets for the dynamics-independent ratio of 
magnetic moments and pion couplings |P[ 

/,AA ' ' -0.71 vs. 0.70 ±0.01 (exp.), 



^^ = ^ = 1.5 vs. 1.5 ±0.12 (exp.). (4.48) 

9-kNN 2 

I should mention that there might be interesting implications of the 'baryons as rotating 
solitons' idea to nuclear physics. The low-energy interactions between nucleons can be viewed 
as interactions between spherical tops depending on their relative orientation -R1-R2 in the 
spin-isospin spaces |55|, [56|1 . It leads to an elegant description of NN and NA interactions 
in a unified fashion, and it would be very interesting to check its experimantal consequences 
(as far as I know this has not been done yet). A nuclear medium is then a medium of 
interacting quantum spherical tops with extremely anizotropic interactions depending on 
relative orientations of the tops both in the spin-isospin and in ordinary spaces. 

This unconvential point of view is strongly supported by the observation [|5l| that one can 
get the correct value of the so-called symmetry energy of the nucleus, 25 MeV ■ (N — Z) 2 /A, 
with the coefficient 25 MeV appearing as 1/ 8/1 where I\ is the SU (2) moment of inertia; 
from the A — iV splitting Q4.47 ) one finds I\ ~ (200 MeV)^ 1 . I do not know whether the 



language of spherical tops is fruitful to describe ordinary nuclear matter (probably it is but 
nobody tried), however it is certainly useful to address new questions, for example whether 
nuclear matter at high densities can be in a strongly correlated antiferromagnet-type phase 

II- 

Finally, let us ask what the next rotational excitations could be? If one restricts oneself 
to only two flavours, the next state should be a (5/2, 5/2) resonance; in the three-flavour 
case the third rotational excitation is the anti-decuplet with spin 1/2, see above. Why do not 
we have any clear signal of the exotic (5/2, 5/2) resonance? The reason is that the angular 
momentum J = 5/2 is numerically comparable to N c = 3. Rotations with J « N c cannot 
be considered as slow: the centrifugal forces deform considerably the spherically-symmetric 
profile of the soliton field |T], ^7|; simultaneously at J ~ N c the radiation of pions by the 



rotating body makes the total width of the state comparable to its mass |], [37], |58 |. In order 



to survive strong pion radiation the rotating chiral solitons with J > N c have to stretch into 
cigar-like objects; such states lie on linear Regge trajectories with the slope a' ~ 1/8ti 2 F 2 
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The situation, however, might be somewhat different in the t/iree-flavour case. First, the 
rotation is, roughly speaking, distributed among more axes in flavour space, hence individual 
angular velocities are not neccessarily as large as when we consider the two-flavour case with 
J = 5/2. Actually, the SU(2) baryons with J = 5/2 belong to a very high multiplet from 
the SU(3) point of view. Second, the radiation by the soliton includes now K and r\ mesons 
which are substantially heavier than pions, and hence such radiation is suppressed. Actually, 



the anti-decuplet seems to have moderate widths |54| and it is worthwhile searching for the 
predicted exotic states. 

4.5 Some applications 

There exists by now a rather vast literature studying baryon observables in the Chiral Quark- 
Soliton Model []. Baryon formfactors (electric, magnetic and axial), mass splittins, the 
nucleon sigma term, magnetic moments, weak decay constants, tensor charges and many 
other characteristics of nucleons and hyperons have been calculated in the model. I address 
the reader to an extensive review on these matters. 



Here I would like to point out several developments of the Chiral Quark-Soliton Model 
interesting from the theoretical point of view. The list below is, of course, very subjective. 

The study of the spin content of the nucleon in the model has been pioneered by Waka- 
matsu and Yoshiki They showed that the fraction of the nucleon spin carried by the 
spin of quarks is about 50% (and could be made less): the rest is carried by the interquark 
orbital moment, the Dirac sea contribution to it being quite essential. 

An important question is 1/N C corrections to baryon observables. These can be classified 
in two groups: one comes, for example, from meson loops and is therefore accompanied by 
an additional small factor ~ l/8n 2 , the second arises from a more accurate account for the 
quantization of the zero rotational modes. The second-type corrections are not accompanied 
by small loop factors, and may be quite substantial: after all in the real world N c = 3 so a 
30% correction is not so small. Such corrections for certain quantities have been fished out 
in refs. |62 , |63| in the two-flavour case and in ref. |j4| for three flavours. These corrections 



8 Sometimes the model has been called 'the solitonic sector of the Nambu-Jona-Lasinio model'. I have 
many objections to this title. First, Vaks and Larkin have suggested independently and at the same time 
a 4-fermion model to illustrate how symmetry can be dynamically broken in field theory. Therefore, in any 
case I would call it the VL/NJL model. Second, both VL and NJL [[l8) were talking about nucleons 
as fundamental fermions (there were no quarks in 1961), and this is rather far from what we consider now. 
Third, as discussed in section 2, the instanton-induced interactions in contrast to the ad hoc 4-fermion 
interactions correctly reproduce the symmetries: Ua(1) is explicitly broken, while in the iV c = 2 case they 
possess a more wide SU(A) x t/(4) symmetry at Nf > 2 they necessarily are 2iV/-fermion interactions 
and not at all 4-fermion. Fourth, instanton-induced interactions provide a natural UV cutoff, as given by the 
momentum-dependence of the constituent quark mass. For the success of the Chiral Quark-Soliton Model 
it is extremely important that this UV cutoff is much larger than the constituent quark mass itself (actually 
squared), meaning that the size of the constituent quark is much less than the size of the nucleon. In an ad 
hoc 4-fermion model having no obvious relation to QCD one has to impose the UV cutoff by hands. Fifth, 
in an arbitrarily introduced 4-quark interaction model there are no a priori reasons to freeze out all degrees 
of frcdom except the chiral ones. Meanwhile, if one includes the a field into the minimization of the nucleon 
mass the soliton collapses |S9| pQ] . In short, when one knows results coming from instantons, it is possible 
to mimic some of them by imposing certain rules of the game with the 4-quark interactions. But why then 
should it be called the 'NJL model' ? 
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work in a welcome direction: they lower the fraction of nucleon spin carried by quark spins 
and increase the flavour non-singlet axial constants. 

A recent development of the model deals with the parton distributions in nucleon. This 
topic deserves however a special subsection. 



4.6 Nucleon structure functions 

The distribution of quarks, antiquarks and gluons, as measured in deep inelastic scattering 
of leptons, provides us probably with the largest portion of quantitative information about 
strong interactions. Until recently only the evolution of the structure functions from a high 
value of the momentum transfer Q 2 to even higher values has been succesfully compared 
with the data. This is the field of perturbative QCD, and its success has been, historically, 
essential in establishing the validity of the QCD itself. However, the initial conditions for 
this evolution, namely the leading- twist distributions at a relatively low normalizaion point, 
belong to the field of non-perturbative QCD. If we want to understand the vast amount of 
data on unpolarized and polarized structure functions we have to go into non-perturbative 
physics. 

The Chiral Quark-Soliton Model presents a non-perturbative approach to the nucleons, 
and it is worthwhile looking into the parton distributions it predicts. Contrary to several 
models of nucleons on the market today, it is a relativistic field-theoretical model. This cir- 
cumstance is of crucial importance when one deals with parton distributions. It is only with a 
relativistic field-theoretical model one can preserve general properties of parton distributions 
such as 



• relativistic invariance, 

• positivity of parton distributions, 

• partonic sum rules which hold in full QCD. 



There are two seemingly different ways to define parton distributions. The first, which I 
would call the Fritsch-Gell-Mann definition, is a nucleon matrix element of quark bilinears 
with a light-cone separation between the quark if) and ip operators. According to the second, 
which I would call the Feynman-Bjorken definition, parton distributions are given by the 
number of partons carrying a fraction x (the Bjorken variable) of the nucleon momentum in 
the nucleon infinite-momentum frame. See Feynman's book |65] for the discussion of both 
definitions. In perturbative QCD only the Fritsch-Gell-Mann definition has been exploited 
as one has no idea how to write down the nucleon wave function in the infinite-momentum 
frame, which is necessary for the Feynman-Bjorken definition. 

Despite the apparent difference in wording, it has been shown for the first time, whithin 
the field-theoretical Chiral Quark-Soliton Model, that the two definitions are, in fact, equiv- 
alent and lead to identical working formulae for computing parton distributions: in ref. |66 



the first definition has been adopted while in ref. |)7| the second was used. The deep reason 
for that equivalence is that the main hypothesis of the Feynman-Bjorken parton model, 
namely that partons transverse momenta do not grow with Q 2 [^5| , is satisfied in the model. 
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Let me point out some key findings of refs. pBI, B7 



(i) Classification of quark distributions in N c 



Since the nucleon mass is 0(N C ) all parton distributions are actually functions of xN c . 
Combining this fact with the known large- N c behaviour of the integrals of the distributions 
over x one infers that all distributions can be divided in 'large' and 'small'. The 'large' 
distributions are, for example, the unpolarized singlet and polarized isovector distributions, 
which are of the form 

D 1 ^ c (x)^N 2 J(xN c ), (4.49) 

where f(y) is a stable function in the large- N c limit. On the contrary, the polarized singlet 
and unpolarized isovector distributions give an example of 'small' distributions, having the 
form 

D^ix) ~ N c f{xN c ). (4.50) 

One, indeed, observes in experiment that 'large' distributions are substantially larger than 
the 'small' ones. 



(ii) Antiquark distributions 



In the academic limit of a very weak mean pion field in the nucleon the Dirac continuum 
reduces to the free one (and should be subtracted to zero) while the bound-state level joins 
the upper Dirac continuum. In such a limit there are no antiquarks, while the distribution 
of quarks becomes q(x) = N 2 S(xN c — 1). In reality there is a non-trivial mean pion field 
which a) creates a bound-state level, b) distorts the negative-energy Dirac continuum. As the 
result, the above 5-functiuon is smeared significantly, and a non-zero antiquark distribution 
appears. 

An inevitable consequence of the relativistic invariance is that the bound-state level 
makes a negative-definite contribution to the antiquark distribution |. The antiquark distri- 
bution becomes positive only when one includes the contribution of the Dirac continuum. 
Numerically, the antiquark distribution appears to be sizeable even at a low normalization 
point, in accordance with phenomenology. 



(iii) Sum rules 



The general sum rules holding in full QCD are automatically satisfied in the Chiral 



Quark-Soliton Model: in refs. |66j, 67] the validity of the baryon number, isospin, total 



momentum and Bjorken sum rules has been checked. In fact, it is for the first time that 



9 This is also true for any nucleon model with valence quarks, for example for any variant of bag models. 
Bag models are essentially non-relativistic, so they fail to resolve this paradox. In order to cure it, one has 
to take into account contributions to parton distributions from all degrees of freedom involved in binding 
the quarks in the nucleon. That can be consistently done only in a relativistic field-theoretical model, like 
the one under consideration. 
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nucleoli parton distributions at a low normalization point have been consistently calculated 
in a relativistic model preserving all general properties. 



(iv) Smallness of the gluon distribution 



As many times stressed in these lectures, the whole approach of the Chiral Quark-Soliton 
Model is based on the smallness of the algebraic parameter (Mp) 2 where p is the average 
size of instantons in the vacuum. This p is the size of the constituent quark, while the size 
of the nucleon is, parametrically, 1/M. Computing parton distribution in the model one is 
restricted to momenta k <C 1/p « 600 MeV, so that the internal structure of the constituent 
quarks remains unresolved. There are no gluons in the nucleon at this resolution scale; 
indeed, the momentum sum rule is satisfied with quarks and antiquarks only. 

However, when one moves to the resolution scale of 600 MeV or higher, the constituent 
quarks cease to be point-like, and that is at this scale that a non-zero gluon distribution 
emerges. Having a microscopic theory of how quarks get their dynamical masses one can 
compute the non-perturbative gluon distribution in the constituent quarks [^ What can be 
said on general grounds is that the fraction of momentum carried by gluons is of the order of 
(Mp) 2 « 1/3, which seems to be the correct portion of gluons at a low normalization point 
of about 600 MeV where the normal perturbative evolution sets in. 



(v) Comparison with phenomenology 



There are several parametrizations of the nucleon parton distributions at a relatively low 
normalization point, which, after their perturbative evolution to higher momentum transfer 
Q 2 , fit well the numerous data on deep inelastic scattering. The most daring (and convenient 
for our purpose) parametrization is that of Gluck, Reya et al. [^, ^] who pushed the nor- 
malization point for their distributions to as low as 600 MeV, starting from the perturbative 
side. In refs. pB] , |6T| parton distributions following from the Chiral Quark-Soliton Model 
have been compared with those of refs. |69|, f70"| . There seems to be a good qualitative agree- 
ment though the constituent quark and antiquark distributions appear to be systematically 



'harder' than those of |T0|| . This deviation is to be expected since the structure of the 
constituent quarks themselves has not been yet taken into account, see above. 



5 Conclusions 

The Chiral Quark-Soliton Model is a simple and elegant reduction of the full-scale QCD at 
low energies, however preserving its main ingredients, namely spontaneous chiral symmetry 
breaking, and the appearance of the dynamical (or constituent) quark mass. Personally, 
I prefer the word 'dynamical': first, because it is, indeed, dynamiclly generated, second, 
because it is momentum-dependent. 

The momentum dependence of the dynamical quark mass M(k) is the key to understand- 
ing why the notion of constituent quarks have worked so remarkably well over 30 years in 
hadron physics. The point is, the scale A at which the function M(k) falls off appears to 

10 Steps in that direction has been taken in refs. |H| |68j. 
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be much larger than M(0); the former parameter determines the size of constituent quarks 
while the latter parameter determines the size of hadrons. These two distinctive scales come 
neatly from instantons, as described in section 2. 

The Chiral Quark-Soliton Model fully exploits the existence of the two distinctive scales: 
it is because of them it makes sense to restrict oneself to just two degrees of freedom in the 
nucleon problem, namely to massless or nearly massless (pseudo) Goldstone pions and to 
the constituent quarks with a momentum-dependent dynamical mass M(k). The scale A 
actually plays the role of the physical ultra-violet cutoff for the low-energy theory; its domain 
of applicability is thus limited to the range of momenta k ~ M < A. This is precisely the 
domain of interest for the nucleon binding problem. 

A technical tool simplifying considerably the nucleon problem is the use of the large N c 
logic. At large N c the nucleon is heavy, and one can speak of the classical self-consistent 
pion field binding the N c valence quarks of the nucleon together. The classical pion field 
(the soliton) is found from minimizing the energy of the bound-state level plus the aggregate 
energy of the lower Dirac continuum in a trial pion field. The valence quarks (sitting on the 
bound-state level) appear to be strongly bound by the classical pion field. 

By quantizing the slow rotations of the soliton field in flavour and ordinary spaces one 
gets baryon states which are rotational excitations of the static 'classical nucleon'. The 
classification of the rotational excitations depends on the symmetry properties of the soliton 
field, but not on the details of dynamics. Taking the hedgehog ansatz one gets the following 
lowest baryon multiplets: octet with spin 1/2, decuplet with spin 3/2 (these are, indeed, the 
lowest multiplets observed in nature) and antidecuplet with spin 1/2. This last multiplet 
contain baryons with exotic quantum numbers (in the sense that they cannot be composed of 
only three quarks); some of them are predicted to be relatively light and narrow resonances, 
and it would be of great interest to search for such states. 

By sayng that all lightest baryons are nothing but rotational excitations of the same 
object, the 'classical nucleon', we get many relations between members of baryon multiplets 
which are all realized with astonishing accuracy in nature. Especially successful are predic- 
tions which do not depend on dynamical quantities (like the values of moments of inertia) 
but follow from symmetry considerations only, and are therefore shared, e.g., by the Skyrme 
model. Predictions of the Chiral Quark-Soliton Model which do depend on concrete dynam- 
ics are, in general, also in good accordance with reality: the typical accuracy for numerous 
baryon observables computed in the model is about 15-20%, coinciding with the expected 
theoretical accuracy of the model. To get a better accuracy one needs a better understanding 
of the underlying QCD vacuum and of the resulting effective low-energy theory. The devel- 
oped theory of the instanton vacuum of QCD seems to do the job of explaining the hadron 
world pretty well already, however if one wants to improve the accuracy of predictions one 
has to make the theory more precise. 

To my knowledge, the Chiral Quark-Soliton Model is the only relativistic field-theoretical 
model of the nucleon on the market today, and this advantage of the model becomes crucial 
when one turns to the numerous parton distribution functions. It is impossible to get a 
consistent description of parton distributions satisfying positivity and sum rules restrictions, 
without having a relativistic theory at hand and without taking into account the complete set 
of forces which bind quarks together. The leading-twist parton distributions computed so far 
in the Chiral Quark-Soliton Model refer to a very low normalization point where the structure 
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of the constituent quarks is not resolved yet. Nevertheless, they seem to be in qualitative 
agreement with parametrizations of the DIS data at low Q 2 though, not unnaturally, they 
appear to be more 'hard'. 

I think that it is the field of parton distributions where the Chiral Quark-Soliton Model 
will be used most of all in the near future. We know how to (perturbatively) evolve parton 
distributions fron high to still higher values of Q 2 but we do not really understand how to 
explain the initial conditions for that evolution, that is the leading-twist parton distribu- 
tions at a low normalization point. This is where a relativistic model satisfying all general 
requirements could be of great use. Also in the years to come there will be much exper- 
imental activity involving numerous spin and off-forward parton distributions, as well as 
non-leading-twist distributions. Practically nothing is known about these numerous distri- 
butions from the theoretical side, and the predictions of the Chiral Quark-Soliton Model can 
be very valuable, see refs. [BR W\ [72[ for the first predictions. 
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